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Abstract 



We employ the QCD sum rules method for description of micleons in nuclear matter. 
We show that this approach provides a consistent formalism for solving various problems of 
nuclear physics. Such nucleon characteristics as the Dirac effective mass m* and the vector 
self-energy Sy are expressed in terms of the in-medium values of QCD condensates. The 
values of these parameters at saturation density and the dependence on the baryon density 
and on the neutron-to-proton density ratio is in agreement with the results, obtained by 
conventional nuclear physics method. The contributions to m* and are related to 
observables and do not require phenomenological parameters. The scalar interaction is 
shown to be determined by the pion-nucleon cr-term. The nonlinear behavior of the scalar 
condensate may appear to provide a possible mechanism of the saturation. The approach 
provided reasonable results for renormalization of the axial coupling constant, for the 
contribution of the strong interactions to the neutron-proton mass difference and for the 
behavior of the structure functions of the in-medium nucleon. The approach enables to 
solve the problems which are difficult or unaccessible for conventional nuclear physics 
methods. The method provides guide-lines for building the nuclear forces. The three- 
body interactions emerge within the method in a natural way. There rigorous calculation 
will be possible in framework of self-consistent calculation in nuclear matter of the scalar 
condensate and of the nucleon effective mass m*. 

PACS numbers 21.65.-f, 21.65.Mn, 24.85. +p 
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1 Introduction 

In the present paper we review our approach to description of a nucleon, placed into nuclear 
matter. The in-medium characteristics of nucleon can tell us much about the medium itself. 
On the other hand, the introduction of nuclear matter enables to separate the problems of 
nucleon interactions from those, connected with individual features of specific nuclei. Thus 
investigation of nucleon characteristics in nuclear matter is an important step for studies of 
physics of nuclei. 

Until mid 70-th the studies of nuclear matter were based on the nonrelativistic approach. 
Since the publication of the paper [1] description of the in-medium nucleon was based on Dirac 
phenomenology. It was successful in describing most of characteristics of nucleons both in 
nuclear matter and in finite nuclei [21 E] • In the meson exchange picture the vector and scalar 
fields correspond to exchange by the vector and scalar mesons between a nucleon and the 
nucleon of the matter. This picture is called Quantum Hadrodynamics (QHD). In the simplest 
version (QHD-I) only vector u and scalar a mesons are involved. In more complicated versions 
some other mesons are included [21 |4j. While the studies in framework of the nonrelativistic 
approach are going on, i.e. the applications of the Nuclear Density Functional Method provided 
very accurate description of the data [3 [6]. 

On the other hand, many efforts have been made to improve the QHD. There were many 
reasons for this. First of all, the QHD itself has several weak points. It is not clear, if the 
scalar a meson does exist, the experimental data are controversial. It is rather an effective 
way of describing the two-pion exchange. The mass of this effective state is about 500 Mev. 
The mass of vector u meson is about 780 MeV. Hence, the exchange by theses mesons takes 
place at the distances, where the nucleon can not be treated as a point -like particle. Thus the 
QHD inherited the problems of nonrelativistic phenomenology, connected with description of 
the interaction at small distances. The other weak points are discussed in [71 E]. Also, it is 
desirable to match the QHD with the Quantum Chromodynamics (QCD ) which is believed to 
be a true theory of strong interactions. Another reason is that there are various problems in 
nuclear physics. It is desirable to have an approach, which would enable to calculate: 

• The nucleon single-particle potential energy f/(p), where p in the density of the baryon 
quantum number. This enables one to find the saturation density po and the single- 
particle binding energy e{p). 

• Parameters of interaction with external fields. These are magnetic moments /i(p) and 
the axial coupling constant 5'a(p)- The latter is important for understanding of the chiral 
properties of the matter. 



• Neutron-proton mass splitting in isotope-symmetric matter. It was observed for a number 
of nuclei, being known as the Nolen-Scliiffer anomaly. 

• Structure functions of the deep inelastic scattering. They describe the internal structure 
of nucleons. Investigation of the latter is important for construction of the quark models 
of nucleons, for studies of the confinement, etc. The reasons for some of the in- medium 
modifications of the structure functions are still obscure. They become a subject for 
discussions from time to time. 

• The single particle potential energy for hyperons in nuclear matter. Can a system of 
hyperons be stable? In other words, can a strange matter exist ? 

As we said earlier, the first problem was solved in framework of Dirac phenomenology 
[H [2] for the values of density close to the saturation value. The nucleon was considered as a 
relativistic particle, moving in superposition of vector and scalar fields and $. In the rest 
frame of the matter = Vo6f^o. The dynamics of the nucleon is described by equation 

{q-V)tlJ = {m + ^)ilj, (1) 

with A = A^7^, = —id^. In nuclear matter the fields V and $ depend only on the density p, 
and do not depend on the space coordinates. The values of the fields are adjusted to reproduce 
either the data on nucleon-nucleon scattering or the nuclear data - see, e.g. |3]. However, each 
of the other listed problems requires additional improvements of the QHD. Turning to the other 
problems from the above list, note that the axial coupling constant changes due to polarization 
of medium by pions [9], with a crucial role of the delta-isobar excitations [10]. Thus, in order 
to solve the second problem one should introduce additional degrees of freedom into the QHD. 
Also, the third problem requires introduction of the nuclear forces, which break the isospin 
invariance [TT]. Finally, the fourth problem is just unaccessible for traditional methods of 
nuclear physics. 

There were several attempts to combine the QHD and the quark structure of nucleons. The 
nucleon-nucleon forces, constructed within such models enabled to reproduce semi-quantitatively 
or quantitatively the nucleon characteristics in nuclear matter. In the Quark-Meson Coupling 
(QMC) model [12] the nucleon was considered as a three-quark system in a bag. The quarks 
were directly coupled to a and u mesons. The values of the nucleon effective mass m* and of 
the in-medium coupling constant qa appeared to be somewhat smaller than in QHD |13j. This 
model is employed nowadays as well |14j . 

In another class of works the short-range interaction between the nucleons was treated 
as interaction between their quarks. The latter were described in framework of some QCD 
motivated models. The long-range NN interaction was described in terms of nucleons and 
pions. These works were reviewed in |15j . 

During the two latest decades much work was done on development of the Effective Field 
Theory (EFT). The starting point is the most general Lagrangian, which includes nucleons 
and pions as the degrees of freedom and respects all the symmetries of QCD, i.e. the Lorentz 
invariance, chiral symmetry, etc. [12]. The applications of the EFT is usually combined with 
the expansion in powers of the pion mass (Chiral Perturbation Theory). Expansion in powers 
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of low momenta is also usually carried out. The works, based on the EFT are reviewed in |17] . 
for some of recent developments - see [IB]. The EFT approach does not touch the quark degrees 
of freedom of the nucleons. 

All the traditional nuclear physics approaches face difficulties in attempts to describe the 
nucleon-nucleon (NN) interaction at small distances. On the other hand, the Quantum Chromo- 
dynamics (QCD) is believed to be a true theory of strong interactions. It has many unsolved 
problems at large distances. However, it becomes increasingly simple at small distances due 
to the asymptotic freedom. The interaction of quarks and gluons, which are the ingredients 
of QCD can be treated perturbatively. This is known as the asymptotic freedom [19]. It is 
tempting to use this feature of QCD in building the nucleon forces. One should take into 
account, however, that due to spontaneous breakdown of the chiral symmetry of the QCD, the 
vacuum expectations of some QCD operators (condensates) have nonzero values. 

In medium with a nonzero value of the density of the baryon quantum number the con- 
densates change their values. Also, some other condensates, which vanish in vacuum, obtain 
nonzero values. 

The QCD sum rules (SR) method describes the vacuum hadron parameters basing on the 
quark dynamics at short distances, where the asymptotic freedom works. In other words, the 
dynamics of the quark system at the distances of the order of the confinement radius is described 
basing on that at small distances, where it is determined by the QCD condensates. Thus, the 
SR method enables to describe the hadron parameters in terms of the QCD condensates. 

The approach was worked out in [20], where it was used for mesons. It is based on the 
dispersion relation for the function, describing the system which carries the quantum numbers 
of the hadron. The SR method was successfully applied for nucleons in vacuum [21], describing 
all their static and some of dynamical characteristics - see [221 123] for a review. Thus it looks 
reasonable to try to apply the QCD SR method for the description of nucleons in nuclear 
matter. It was suggested in [211 |25l [26] that the parameters of nucleon in nuclear matter can 
be expressed in terms of the in-medium values of QCD condensates. In medium with a nonzero 
value of the density of the baryon quantum number the condensates change their values. Also, 
some other condensates, which vanish in vacuum, obtain nonzero values. 

The generalization of the SR method for the case of finite densities was not straightforward. 
One of the main problems was the choice of variables, which enabled to separate the singularities 
connected with the in-medium nucleon from those connected with the medium itself. This was 
done in [21]- [Si- 
lt was found also in these papers that the nucleon characteristics (effective mass m* and 
the vector self-energy Ey) can be presented in terms of the vector and scalar condensates 



Here p and (M| are the density and vector of the ground state of the matter, is the quark 
field, the summation over flavors i = u,d is carried out. The vector condensate is written in 
the rest frame of the matter. Due to conservation of the vector current the vector condensate 
is a linear function of p 



v{p) = {M\J2<fjoq'\M)- K{p) = {M\Y,<fq^M). 



(2) 



v{p) = VnP; = {N\ ^ q'joq'lN) = 3 



(3) 



is just the number of valence quarks in a nucleon. The scalar condensate can be represented as 

[MIES] 

K{p) = k{Q) + KMP + 5(p), KM = {N\ q'q'\N) , (4) 

i 

with S{p) caused by interaction of the nucleons of the matter. Since can be expressed 
in terms of the pion-nucleon a term [27], while the latter is related to observables [281 129] . 
one can obtain the values of the nucleon parameters at least in the gas approximation. Since 
S[p) is small for the densities, close to the saturation point (see below), the values obtained 
in such approach are close to the physical ones. They appeared to be Sy ~ 200 MeV, and 
m* —m^ —300 MeV close to the saturation point. Thus, including only the condensates of the 
lowest dimension and neglecting the radiative corrections, we found that the QCD SR method 
reproduce the main features of the QHD [T]. 

However, the role of the condensates of higher dimension remained obscure. The contribu- 
tion of the gluon condensate (Ml^G^'^^GjJj^lM) is rather small. However, an estimation for the 
value of the four-quark condensate {M\qqqq\M) , which suggests itself, destroys the agreement 
with the Walecka model. Also, the lowest order radiative corrections are numerically large. 
This took place for the vacuum as well. This caused doubts in possibility to expand the SR 
method for the case of finite densities. 

The four-quark condensates are the most important among those of the higher dimension. 
Their calculation requires some model assumption on the quark structure of nucleon. We 
employed the Perturbative Chiral Quark Model (PCQM), suggested originally in [30]. We 
calculated these condensates |31j and found that previous naive estimation of its contribution 
was wrong, due to some cancelations which take place in any reasonable model. We demonstrated 
that the SR method provides the dependence of the nucleon characteristics m* and on the 
density p and on the neutron-to-proton density ratio [32l [33], which is consistent with the 
results obtained by using traditional nuclear physics methods. 

We analyzed the role of radiative corrections for the nucleon SR in vacuum and demonstrated 
that their influence on the value of the nucleon mass is small [34]. Also, in nuclear matter the 
radiative corrections do not change much the nucleon characteristics m* and Ey [HS]. 

We found that the nonlinear contribution to the scalar condensate S{p) is determined mainly 
by the pion contribution to the self-energy of the nucleon of the matter. Simple estimations show 
that this term may provide a saturation mechanism in our approach [36]. However, a rigorous 
treatment requires renormalization of the pion propagator by the particle-hole excitations. The 
renormalized pion propagator depends on the nucleon effective mass m* and on the in-medium 
value of the pion-nucleon couphng constant ^f^.The latter can be obtained by the SR method. 
This brings us to a self-consistent scenario. At the present level of our knowledge it requires 
some more phenomenological assumptions [37] . 

Note also, that the finite density SR method provided reasonable results for the axial cou- 
pling constant [3H], for the neutron-proton mass splitting [5U] and for the difference between 
the deep inelastic structure of nucleus and that of sum of those of free nucleons [IQ] . However, 
in these calculations only the condensates of the lowest dimensions have been included. 

Now we give the details. 



2 Nucleon QCD sum rules in vacuum 



This approach is described in details in many pubhcations. Earher papers are reviewed in ^T\ . 
The nowadays state of art is presented in [22] and [23], see also the book [19]. However, to 
make the text self-consistent we recall the main points of the approach. We emphasize the 
points, which we shall need for the extension of the SR method for the case of finite density of 
the baryon quantum number. 

2.1 General ideas 

The nucleon QCD sum rules succeeded in describing of the nucleon characteristics in vacuum 
in terms of the vacuum expectation values of the products of quark or (and) gluon operators 
(QCD condensates) [3 [8]. This approach is based on the dispersion relation for the function 



describing the propagation of the system which carries the quantum numbers of the proton, / 
is the unit 4x4 matrix. In the simplest form the dispersion relations are 



As we shall see below, we do not need to worry about possible subtractions. 

In quantum mechanics no(g^) is just the proton propagator. In the field theory different 
degrees of freedom are important in different regions of the value of g^. 

One can consider the proton as a system of three strongly interacting quarks. Due to 
asymptotic freedom of QCD 

citelSa the description becomes increasingly simple at — ?■ — oo. This means that at — t- 
— oo the function no(g^) can be presented as a power series of and of the QCD coupling 
constant as- The coefficients of the expansion in powers of are the QCD condensates. Such 
presentation known as the Operator Product Expansion (OPE) [12] provides the perturbative 
expansion of the short distance effects, while the nonperturbative physics is contained in the 
condensates. In QCD SR approach the left-hand side (LHS) of Eq.(l5]) is considered at — )■ — oo, 
and several lowest order terms of the OPE are included. In this and in the next Subsections 
we neglect the radiative corrections, i.e. we do not include interactions and self-interactions 
(self-energy insertions) of the quarks, putting = 0. 

Turning to the right-hand side (RHS) of Eq. ([5]) note that Imno(A;^) = at /c^ < 
with m being the position of the lowest lying pole, i.e. m is the proton mass. There are 
the other singularities at larger values of fc^. These are the cuts corresponding to the systems 
"proton-|-pions" , the pole N(1440), etc. The next to leading singularity is the physical branching 
point k"^ = Wpi^ , and one can write 




(5) 



with = 1, e 




(6) 
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k'^ — q'' 



dk 
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(7) 
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The upper index OPE means that several lowest OPE terms are included. 

The SR approach is focused on studies of the lowest state. Thus we keep the first term on 
the RHS of Eq. ([7j) and try to write approximate expression for the contribution of the higher 
states. The detailed structure of the spectral function /(fc^) in the second term on the RHS of 
Eq. ([7]) cannot be obtained by the SR method. However, at k"^ ^ 

ATTi OPEfu2\ 



2^ 

with A denoting the discontinuity. 

The standard ansatz consists in extrapolation of Eq. (8) to the lower values of k"^, assuming 
that the cut starts at certain unknown point W"^. In other words 



vr _ q2 ^-ni Jw'^ k"^ — 

vnys 



and thus Eq. (5) takes the form 



\2 H 1 ATT* OPEtv.'i^. 



rn? — 2m J k"^ — q"^ 

Recall that ^-J = 1, = m. 

This "pole+continuum" presentation of the RHS makes sense only if its first term, treated 
exactly is larger than the second term, which approximates the higher states. The position of 
the lowest pole m, its residue and the continuum threshold are the unknowns in Eqs. 

In the next step one usually applies the Borel transform. It is defined as 



BFm = Jim =^^(-3;^ W)-F(M^ 



Q' = -q', M' = Q'/n, (10) 

converting a function of into the Borel transformed function of the Borel mass M^. The 
reasons for applying the transform are 

• Since B[Q'^Y = ^"^Y integer k, it kills the polynomials of q^. Hence, it eliminates the 
divergent terms in the function Hq. This explains, why we wrote the dispersion relation 
(5) without subtractions. 

• It emphasizes the contribution of the lowest state, since 

Q2 _j_ ^2 

• It improves the convergence of the OPE series, since 
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The Borel transformed SR take the form 

III '''''' {M^) = A^f e(-"''/^') + — / de^II^ «^^(A;2)e(-^'/^^'). (11) 

2711 J 

If both LHS and RHS of Eq. fllip were calculated exactly, the relation would be independent 
on M^. However, certain approximations are made on both sides. The analytical dependence 
of both sides on is quite different. The OPE on the LHS becomes increasingly true at large 
M^. The accuracy of the "pole+ continuum" model used for the RHS increases at small M^. 
An important assumption is that there is an interval of the values of M^, where the two sides 
have a good overlap, approximating also the true functions Hq(M^). 

Thus our task is to find the region of the values of M^, where the overlap can be achieved 
and to obtain the characteristics of the lowest state m and and the value of the continuum 
threshold W"^. 



2.2 Explicit form of the SR equations 

The general definition of the function Ho(g^), which is sometimes called "correlator" or the 
"correlation function" is 

Ho(g2) = I I d^a;e*(«-^)(0|T[j(x)j(0)]|0), (12) 



where j is the local operator with the proton quantum numbers. It is often refereed to as 
"current". It was shown in |36] that there are three independent currents 

ji = {ulC^^u,)^,rdce'''"^ 32 = {ulCcj^,.Uhh5a''^''d,e^''', (13) 

Here a, b, c are the color indices, C is the charge conjugation matrix, while T denotes the 
transpose in the Dirac space. It was shown in ^3] that the operators j2 and ja^ provide strong 
admixtures of the states with negative parity and of the states with spin 3/2 correspondingly. 
Thus, the calculations with the operator ji are most convincing. We shall assume j = ji in 
further studies. 

Expansion of the matrix element on the RHS of Eq. (1121) in powers of x'^ corresponds to 
expansion of its LHS in powers of q~^. In the lowest orders of expansion the T product on 
the RHS of Eq. f fT2|) can be written in terms of those of two quark operators. The latter can 
be written, following the Wick theorem 

(0|T [gr(x)g)(0)] |0) = g,{x) - ^ E r^'^"'(0|gTMo) + 0{x'). (14) 



Here 

g,{x) = _l_ (^-^yV2).^ ., ^^5) 



is the free quark propagator, nig is the quark mass, F are the basic 4x4 matrices with the 
scalar, pseudoscalar, vector, pseudovector and tensor structures, i.e. /, 75, 7^, 757^ and aap- For 



P 



the fields which respect the chiral invariance all the expectation values on the RHS of Eq. fll4p 
vanish. However in the QCD the expectation value (0|gg|0) {Ta = I) has a nonzero value. All 
the other condensates (0|gr^g|0) vanish due to invariance of vacuum. Thus 

(0|r [qt{x)<f^{0)] |0) = g,{x) - l/,/'^''(0|gg|0) + 0{x'). (16) 

Now we present 

n n 

with n denoting the dimension of the condensate, contributing to the term An or 5„. 

If all expectation values of the T products of the quark fields are described by the free prop- 
agators (IT^ . we find the leading OPE contribution Aq to the structure no(g^), corresponding to 
the free three-quark loop. If one of the quark pairs is described by the second term on the RHS 
of Eq. (14), while the others are given by the first term, we find the leading contribution to 
the structure HQ(g^). In this case two quarks form a free loop, while the current exchanges by 
a quark-antiquark pair with vacuum. 

Direct calculation provides [2T] 

_ q'H-q^lLl) _ (OM'rflO)g^ln(-gVL^) 

' . (17) 

Here Lq is the cutoff of the integral over x in Eq. (IT^ . Its value is not important, since the 
terms containing Lq will be eliminated by the Borel transform. 

Going beyond the single-particle presentation f|T^ one finds also the next to leading OPE 
terms. For example, the contribution is due to the lowest order interaction of the quark 
system with the gluon condensate 

-n ^ TT 

with E and B the color-electric and color-magnetic fields. This condensate also has a nonzero 
value only due to the violation of the chiral symmetry in the ground state of the QCD. The 
higher terms contain the four-quark condensate 

2(0|M|0) .... 
^6 (18) 

Finally, the Borel-transformed sum rules can be written as 

Cl{M^, W^) = W{M^)- Ci{M\ W^) = R\M^). (19) 

Here 

W{M^) = X\-^"'^'"; R\M^) = mA^e-™'/*^', (20) 

with = 327r^A^. The factor 327r^ is introduced in order to deal with the values of the order 
of unity (in GeV units). The contribution of continuum is moved to the LHS of Eqs. ( IT9l) (see 
Eq. (l22l) below). Following [211 Hi] we can write 

= Aq{M^ ,W^) + Ai^m^ ,W^) + Aq{M^)] = B^{M^ ,W^) + Bj{M^) . (21) 



in 



The terms An and i?„ are the Borel transforms of the contributions An and -B„ to the functions 
Ul'^q"^), with subtraction of the corresponding contributions of continuum. Recall that the 
lower index denotes the dimension of the condensate. The terms, proportional to condensates 
of higher dimension Aq and Bj do not contain the logarithmic loops and thus do not contribute 
to continuum. Actually the calculations are carried out in the chiral limit rriq = 0. The explicit 
form of the contributions is [2T1 HI] 



r 4 ' 3 ' 



^3 = 2aM^E, ; ^7 = • (22) 

Here 

Eo(x) = 1 - e""^, Ei(x) = 1 - (1 + x)e-'^, E2{x) = I - {I + x + /2)e-'' , (23) 

while 

a = -(27r)2(0|gg|0); 6 = (27r)2(0|^G"^^G;jO). (24) 

The contributions are illustrated by Fig. 1. Note that we provided these equations mostly as 
illustration of the main ideas and did not include several numerically not very important terms. 

The term Aq presents the contribution of the four-quark condensates {Q\qV^qqV^q\Q) , which, 
generally speaking, obtain nonzero values for all structures Va- It is evaluated under the 
factorization approximation [6] 



(0|gr^gr^g|0) = -^((0|gg|0)) 
lb 



(Trr^)2 - iTr(ri) 



One can find numerical values of the main QCD condensates presented by Eq. fl2^ . There 
is the well known Gell-Mann-Oakes-Renner relation (GMOR) for the scalar condensate [15] 



2fW 

(0|uu + rfrf|0) = — " . (25) 
ruu + rud 

Here /tt and tjIt^ are the decay constant and the mass of the vr meson, iriu and are the 
current masses of u and d quarks. Its numerical value is {Q\uu + dd\Q) = 2(— 240MeV)^. The 
value of the gluon condensate (0|^G'^'^''G^^|0) ^ (0.33 GeV)'' was extracted from the analysis 
of leptonic decay of p and mesons |16]. This data were supported by the QCD sum rules 
analysis of charmonium spectrum |20] . 

Now one must find the set of parameters m, A^, W"^, which insure the most accurate approx- 
imation of 0{M'^) by the functions R^{M'^) and also the interval of the values of M^, where 
this can take place. The set of parameters m, A^, W"^, which minimize the function 

X{m,\ ,W ) = ^\^ [ ^-—2- , (26) 

j i=q,I V j ) 

will be referred to solution of Eqs. (fTTj) . 
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The appropriate interval 

0.8 GeV^ <M^ < lAGeV^ (27) 
and the values of nucleon parameters 

m = 0.93 GeV; = 1.8 GeV^; = 2.1 GeV^ (28) 

were found in [2n HI]. If we fix m = 0.94 GeV, the SR provide 

= 2.0 GeV\ = 2.2 GeV^. (29) 

These values will be employed in the present paper. 

It was shown also in [2JJ that the nucleon mass vanishes if there is no spontaneous breakdown 
of the chiral symmetry, e.g. if (0|gg|0) = 0. Numerically [2T| HI] 

^ -2(27r)^(0|gg|0) . (30) 

In the QCD SR approach the mass of the nucleon is formed due to exchange by quarks between 
the nucleon and vacuum. 

The role of instantons in QCD SR was analyzed in |17], [H]. It was shown that for the 
current ji the contribution of instantons vanishes for the q structure of the SR. In the scalar 
structure instantons form a factor, which change the LHS of SR by about 15%. Following [22] 
we include the contribution into uncertainties of the numerical value of the vacuum expectation 
value (0|gg|0). However, a rigorous analysis requires investigation of the dependence of the 
contribution. 

Note that the solutions (28), (29) are not stable with respect to modification of the values 
of the condensates [IH]. Even for the small changes the absolute minimum of the RHS of 
Eqs. (19) is provided by another solution, i.e. m = 0.6 GeV, A^ = 0.79 GeV^ = 1.0 GeV^ 
[5U] . We treat this solution as an unphysical one, since the contribution of the continuum 
exceeds more than twice that of the lowest pole. This contradicts the key assumption of the 
"pole+continuum" model for the spectrum - see Eq. (9). 

The nucleon SR with another form of nuclear current were obtained in [51]. In [52] it was 
used also for the description of delta isobars. Further we shall mention some other applications. 



2.3 Inclusion of radiative corrections 

A typical radiative correction is shown in Fig. 2. In the analysis, carried out in [21] HI] the most 
important radiative corrections of the order as In have been included. These contributions 
were summed to all orders of {as In Q^)". This is called the Leading Logarithmic Approximation. 
The LLA corrections are expressed in terms of the factor f53] 



where A ^ 150 MeV is the QCD scale, while /i is the normalization point, the standard choice 
is /i = 500 MeV. 

Al = Ao/L^/'; B; = Bs. (32) 
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Here the upper index r indicates inclusion of the radiative corrections. 

The radiative corrections to the OPE terms of the function Ho have been calculated beyond 
the LLA in the lowest order of the as expansion [M] (see also |55]). The results are 



° 12 TT 2 TT /i2/' 6 H 6 TT 3 7r /i^ / ' ^ ' 



^ 2 71 

The running coupling constant in the one-loop approximation, with inclusion of three light- 
est quarks is 

The calculations in [511 155] have been carried out for as = const. Since the linear momenta in 
the loops, corresponding to radiative corrections are of the order q, it is reasonable to assume 
that in Eq. §M> as = as(Q^), converting to as{M^ ^ IGeV^). Assuming ^^(IGeV^) ^ 0.37 
(somewhat larger values are often used nowadays [56]) we find that the radiative correction to 
the contribution Aq changes its value by about 50%. This "uncomfortably large" correction 
was often claimed as the most weak point of the SR approach [57] . 

In [31] we investigated the influence of radiative corrections on the values of characteristics, 
obtained in framework of the Borel transformed nucleon SR in vacuum. We demonstrated 
that inclusion of the radiative corrections in various ways (the radiative corrections are totaly 
neglected, are included in LLA, are taken into account beyond the LLA in the lowest order) 
alter the value of nucleon mass m by about 5%. The radiative corrections modify mainly the 
value of the nucleon residue. 

Note also that inclusion of the radiative corrections diminishes the role of the unphysical 
solution, mentioned above. Once they are included, minimization of the function defined 
by Eq. (126|1 is provided by the physical solution in a broader interval of the values of the 
condensates Il9l [50] . 



3 QCD sum rules in nuclear matter 
3.1 Choice of the variables 

Now we shall try to use the SR approach for calculation of the nucleon parameters in nuclear 
matter. The propagation of the system which has a four-momentum q and carries the quantum 
numbers of the proton is determined by the equation 

n„ = z / d^xe^(^-^)S(x); E{x) = {M\T[j{x)j{0)]\M), (35) 



with |M) the ground state of the nuclear matter. It is just an analog of Eq. (12). We consider 
nuclear matter as a system of A nucleons with momenta Pi, introducing 

P = (36) 



In the rest frame of the matter P ^ (m, 0). 

The spectrum of the function Ilm{(l, P) is much more comphcated than that of the vacuum 
function no(g^). Our main task is to separate the singularities connected with the nucleon 
in the matter from those connected with the matter itself. Include in the first step only the 
two-nucleon interactions. The singularities connected with the matter manifest themselves as 
singularities in variable s = {P + qY- Thus the separation can be done by considering Ilm{q^, s) 
and keeping s = const. Thus we consider the dispersion relations 

UU,U . II ,3T) 

for the three structures {i = q, P, I) of the function 

n„(g2, s) = qUliq', s) + PC(g^ s) + s). (38) 

Each contribution 11* can be viewed as the sum of the vacuum term IIq and that provided by 
the nucleons of the matter 11^ 

Ul = Ul + IV^; no^ = 0. (39) 

This notation will be used also for the other functions. 
We clarify the value of s putting 

s = 4^2^, (40) 

with Eqf being the relativistic value of the nucleon energy on the Fermi surface. This insures 
that the nucleon pole on the RHS of Eq. fl37|) describes the nucleon, added to the Fermi surface. 
For the analysis, carried out in this section we can neglect the bound, thus putting 

s = 4m^ (41) 

This was the choice of variables in our papers |9]-[TT], [321 [33] and [35]-[l0]. It was used 
also in [58], where the approach was used for calculation of the nucleon-nucleus scattering 
amplitude. 

Note that the vacuum dispersion relation (5) can be viewed as a relativistic generalization 
of the nonrelativistic dispersion relation in time component (energy) go 5 known as the Lehmann 
representation [59]. The latter is based on casuality. It converts into a dispersion relation in 

after being combined with the symmetric relation in negative values of go- The reasoning 
does not work in medium, since the Lorentz invariance is lost. To prove the dispersion relation 
(37) we must be sure of the possibility of the contour integration in the complex g^ plane. A 
strong argument in support of this possibility is the analytical continuation from the region 
of real g^ — ?■ —00. At these values the asymptotic freedom of QCD enables one to find an 
explicit expression for the integrand. The integral over the large circle may have a nonvanishing 
contribution. However, the latter contains only polynomials in g^ which are killed by the Borel 
transform. Thus we consider dispersion relations in g^ to be a reasonable choice. 

On the contrary, dispersion relations in go contain all possible excited states of the matter 
on its RHS. To illustrate the latter point, consider photon propagation in medium (see, i.e. 
[60]). 
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In vacuum the propagator of the photon which carries the energy oj and hnear momentum 
k is Z^o = (i^^ — k'^)~^ . It has a pole at = — A;^ = 0. In medium it takes the form 
Dm = {uj'^e{uj, k) — /c^)~^. Here e{uj, k) is the dielectric function, related to the amplitude of the 
photon scattering on the ingredients of the medium. If the photon energies are small enough, 
dependence of e on can be neglected (this is known as the dipole approximation). However, 
e{(jj, 0) is a complicated function of u. It depends on the eigen energies of the medium. The 
same refers to the function Dm{oj). However, the function Dm{q'^) still has a simple pole at the 
point q"^ = = <^'^{^ — reflecting properties of the in-medium photon. Straightforward 

calculation of the value is a complicated problem. The same is true for the in-medium 
proton. The SR are expected to provide the in-medium position of the pole in some indirect 
way. 

The approach, in which the dispersion relation in go cit fixed value of the three-dimensional 
momentum |q| is the departure point was developed in [61]. It was used for description of 
nucleons [61],|62], delta-isobars [63j and hyperons |6l] in nuclear matter. The results are reviewed 
in [65] . However, the possibility to separate the singularities connected with the nucleon in the 
matter from those connected with the matter itself in this approach looks to us to be obscure. 

Thus we expect the Borel transformed dispersion relations in with fixed value of s to be 
more reliable. 



3.2 Operator product expansion 

Note that the condition s = const enables to use the OPE of the LHS of Eq. (37). Indeed, we 
find 

2(Pg) = s-m^-g^, (42) 
and thus in the rest frame of the matter 

2 2 

s — m — q 
2m 

Hence, 

q"^ /Qo ~^ const = c ~ 2m at |g^| — )■ oo . (43) 
The exponential factor on the RHS of Eq. (l35l) can be written as 

with z - the direction of momentum 1. Hence, the integral in Eq. (37) is 

determined by ~ g~^, ~ c~^, and the function S(x) can be expanded in powers of x^, 
corresponding to expansion of Hm^(g^, s) in powers of g~^. Note that condition (43) is the same 
as that for the validity of the OPE for the structure functions of the deep inelastic scattering 

In the case of vacuum the expansion of quark fields g(x) in powers of x was indeed an 
expansion in powers of x^, corresponding to a power series of g~^ for the vacuum function Ho(g^). 
In medium the fields q{x) can be expanded also in powers of {Px). This leads to expansion in 
powers of {Pq)/q'^ of the function H„i(g^,s), and thus, generally speaking, to infinite number 
of condensates in each OPE term. Fortunately, due to the logarithmic g^ dependence of the 
quark loops, the leading OPE terms contain only finite number of condensates. We shall give 
details in the next Section. 
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3.3 Model of the spectrum 

Now we turn to the RHS of Eq. f l371) . The function defined by Eq. fj35|) can be written as 

= {M\T3{x)m\M) = {MA\3{x)\MA+i){MA+i\m\MA)e{x,) 

- (M^|J(0)|M^_i)(M^_i|j(x)|M^)e(-a;o), (44) 

with \Mk) standing for the system with the baryon number K. Here \M) = \Ma) is the ground 
state, summation over all states with K = A ± 1 is assumed. This equation is illustrated by 
Fig. 3 

The matrix element {Ma+iIJIMa) contains the term (A^|j|0), which adds the nucleon (the 
"probe nucleon") to the Fermi surface of the state \Ma), while the rest A nucleons are spectators. 
If interactions of this nucleon with the matter are neglected, this contribution to H(x) has a 
pole at = m^. Interactions of the probe nucleon with the matter shift the position of the 
pole. In the mean field approximation shown in Fig. 4a the shift does not depend on s. Going 
beyond the mean field approximation we find the Hartree self-energy diagrams (Fig.4b), which 
have singularities in s. Due to condition (41) they do not add new singularities in variable 
to the function Ilm{q'^, s). The exchange (Fock) self-energy diagram is shown in Fig. 4c. 

The matrix element (Ma+iUIMa) contains also the terms {Bi\j\0) with Bi standing for 
the system, containing the nucleon and mesons. The current j{x) creates the nucleon and a 
meson with the mass rrix, which is absorbed by the nucleons of the matter - see Fig. 5. This 
contribution has a cut in complex plane, at g^ > -|- 2mmmx. It has also a cut in s, which 
is not important for us, since s is fixed. 

The matrix element {MA-i\i{x)\MA) contains the terms (i?o|j(0)|A^), with Eq describing a 
system with the baryon quantum number equal to zero. The other A — 1 nucleons are spectators. 
The system Bq can be a set of vr mesons, u meson, etc. These contributions depend on the 
variable u{q^) = {P — qY, providing singularities at m > ml, with the mass of state Bq - see 
Fig. 6. Thus the lowest singularity in u corresponds to the branching point g^ = + 2m^, 
corresponding to the real two-pion state in the u channel. A single-particle meson state with 
the mass rrix generates a pole at g^ = + ?ti^/2. The diagram, shown in Fig. 4c. is also one 
of the contributions, which has singularities in the u channel. 

Note that the antinucleon state corresponding to go = —m generates the pole g^ = 5m^, 
shifted far to the right from the lowest lying state. 

Thus the spectrum of the function Ilm{q'^,s) consists of the pole at g^ = m^, a set of 
higher laying poles, generated by the u channel and a set of branching points. The lowest 
lying branching point is separated from the position of the pole g^ = by a much smaller 
distance than in the case of vacuum (g^ = -|- 2171111^^ in the latter case). Note, however, that 
at the very threshold the contribution is quenched since the vertices contain linear moments of 

intermediate pions. Thus the higher singularities can be considered as separated from the pole 

2 2 

r = mt^- 

The situation becomes more complicated if we include the interaction of the probe nucleon 
with n > 1 nucleons of the matter. The corresponding amplitudes depend on the variables 
Sn = {nP + g)^. This causes the cuts, running to the left from the point g^ = m? . Its 
contribution thus is not quenched by the Borel transform. However, as we shall see in Sec. 8, 
such multinucleon interactions require inclusion of the condensates of high dimensions on the 
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LHS of (37). Such contributions do not have logarithmic loops, thus contributing only to 
the pole terms on the RHS of Eq. (37). Hence, this singularities should be disregarded in 
our approach, and the three-body forces are included in the mean-field approximation in our 
approach. 

To summarize the results of this Section, we use Borel transformed dispersion relations of 
the function Umiq^, s) at fixed s. The OPE is used on the LHS. The "pole-l- continuum" model 
is used for the RHS. Now we have three SR equations 

+ :^r dk'A,.Ul^^^^{e, .)e(-'=V^^); z = q,P, I. (45) 

Here and are the position of the nucleon pole in medium and the value of its residue, 
is the in-medium value of the effective threshold. The meaning of the parameters will 
be clarified in next Section. 



4 Nucleon self-energies in the lowest orders of OPE 

In the Subsections 4.1 — 4.4 we consider the symmetric nuclear matter, with equal number of 
protons and neutrons. Also, due to isotope invariance 

{M\dd\M) = {M\uu\M) = {M\qq\M) . (46) 



4.1 General equations 

Start with description of the nucleon pole. We can write for the inverse nucleon propagator in 
medium 

G],' = {G%)-'-^ (47) 
with G% — {q — m)~^ being the propagator of the free nucleon, while 

E = gE„ + — Ep + Ej (48) 
m 

is the general expression for the self-energy in the nuclear matter. In the kinematics, determined 
by Eq. (41) we find 

r -^m 

with 



For the new position of the nucleon pole we find 

1 + Ey/m 



2 _ {s - vr?)Y.ylm - E^ + m*^ 
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while ^ 

^ = (l-S,)(l + Sv^/m)- ^^^^ 



Thus in Eq. (gSD 

e = l; e^ = -Sy; e' = m*. (53) 



Note that Eqs. (1501) correspond to those for the vector self-energy and for the effective mass 
in nuclear physics - see, e.g. |66]. Also in accordance with definition accepted in nuclear physics 
the scalar self-energy is 

Us = TT^* — fn . (54) 
In the nonrelativistic limit the proton dynamics is determined by the potential energy 

U = Sy + S,. (55) 

Keeping the only terms, which are linear in the density p we find a simple expression for the 
shift of the position of the nucleon pole 

rrim — 171 = 11. (56) 

Instead of Eqs. (19) we have now 

Cl^iM^W^) = R\M'y, t = q,PJ. (57) 

Here £^ are the Borel transforms of the LHS of Eq. (37), multiplied by 32it^ - see the previous 
Section, while 

R'{M^) = m*A^e-™"/*^', (58) 

where is the effective value of the in-medium value of the residue, i.e. = Z ■ 327r''A^^, 
following Eqs. (49) and (52). 

Employing Eqs. (57) and (58) one finds 

4.2 Left-hand sides of the sum rules 

4.2.1 Contribution of the condensates of lowest dimension 

The calculation is based on the presentation of the single-quark propagator in nuclear matter 



(M|T[gf(x)gJ(0)]|M) = ±S^_J^!hp3lL^^^ 

- ^Y.^^''\M\m^\{^)\M) , (60) 
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just analogous to Eq. (14). The bilocal operators in the last term on the RHS of Eq. (60) are 
not gauge invariant. The gauge invariant expression 

q{x) = (l + a^a^a + ^XaX^DaDfi + ...^ q{0), (61) 

with Da standing for the covariant derivatives, provides the infinite series of the local conden- 
sates. Anyway, here we are looking for the condensates of the lowest dimension. Hence, we put 
q{x) = q{0) in the last term of Eq. (60). Thus 



- ^S,jS^'{M\q{0)qmM) - ^^^,S'^'{M\mi,QmM). (62) 

The first term presents the singular part of the propagator. It is just the same as in vacuum- 
Eq.(14). We include it in the chiral limit niq — O. Recall that the calculations in vacuum have 
been carried out in the chiral limit as well. The second term also has the same form as in 
vacuum. However, it includes the operator qq averaged over the ground state of the matter. 
The last term contains a new vector condensate, which vanishes in vacuum. In the rest frame 
of the matter it is 



,1 



1 



V 



M + OpV; < = (iV|g(o)7og(o)|iV) . (63) 



In the leading order of the OPE two of the quarks are described by the free propagators Qq, 
while the rest one is presented by one of the two last terms on the RHS of Eq. (62). Present, 
similar to the vacuum case 

^^(^^«) = E^n; C(?',«) = E^n; C(?^^) = E^n (64) 

n n n 

with n denoting the dimension of the condensates. We find immediately that the contribution to 
can be expressed by Eq.(17) for B3 with the matrix element (0| JdjO) replaced by {M\dd\M). 
Thus 

{M\dd\M)qHn{-qyLl) 
^3 ^2 ■ 

The leading contributions to H'' and H^ are determined by the vector condensate 

Hp = ^ / ^ (i)" + v') + x{x, V- + v') exp (i(g • x))) . 
Direct calculation provides 

where v{p) = 3p is the vector condensate - see Eq. (3). The contributions of these terms to the 
LHS of Eqs. (57) are 
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The functions are defined by Eq. (23). Following our notations 

K^ip) = K{p) - k{0) (68) 

We illustrate Eq. ( 1^ by Fig. 7. It describes the exchange by noninter acting quarks with 
the quarks of the matter. Of course, there are strong interactions between the quarks of the 
condensate. The corresponding contributions to the nucleon pole (the RHS of Eq. fl571) ) account 
for the exchanges by the systems of strongly correlated quarks (mesons) with the same quantum 
numbers. They have standard Lorentz structures, i.e. the terms proportional to v{p) and K,{p) 
contribute to vector and scalar structures of Eq. ([T]) correspondingly. This is illustrated by 
Fig. 8. 

4.2.2 Next to leading OPE terms 

Now we include the contribution of the condensates of dimension d = A. These are the contri- 
butions caused by nonlocality of the vector condensate (M|g(0)7og(x)|M) - see Eq. (60) with 
q{x) defined by Eq. (61) to q and P structures on the LHS of Eqs. (57) and also the contribution 
of the gluon condensate to the q structure. 

Recall that in the case of vacuum there was a contribution A4 caused by the vacuum gluon 
condensate - see Eqs. (22) and (24). In the nuclear matter the vacuum gluon condensate 
should be replaced by its in-medium value. The corresponding contribution is 

Mgp = VT jjj^ gp{p) . (69) 

Here we denoted 

Now we consider the contribution of the nonlocal vector condensate. For the sake of sim- 
plicity we shall present the results in terms of the nucleon matrix elements, assuming thus that 
the matrix element {M\q{0)'-foq{x)\M) is a linear function of p. We can write 

9l{x) = (M|g(0)7og(a;)|M) = ^${((Px), x^) + ^x^m<l>{((Px), x'). (71) 
Presenting [57] 

K,{iP^),^') = £dae-^-^^-^fi,{a,x% (72) 

we can expand j,(a,x^) = ?7afe(a) + 2;^m^^^^(a)/8 -|- 0{x^). Here ri[{a) = //(a,0) is the 
contribution of the valence quark with flavor / to the asymptotics of the nucleon structure 
function rja = t]"^ + 77^, normalized by the condition 



g{p) = m-G'^^'GUM) = giO)+9,ip)- (70) 

IT 



/ daria{a) = 3. (73) 
Jo 
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The moments of the function rjf, can be presented in terms of those of T]a and ^a[2S] , e.g. 



darib{a) = - daaria{a). (74) 
4 Jo 

In the equations, displayed in this Subsection we omit some terms, which are not important 
numerically. 

Direct calculation provides the contribution of the nonlocal condensate to 
Pap = Pil^ + Pif; Pip' = \n^J da[{5a{P ■ q') + 2a'm^)ria{a) + Qm^r^,{a)]p- 



p 



^4? = 7r4- /' daln^li-aiPq') + 2q")r]aia) - 9m^ar]tia)]p, (75) 
biT'^m Jo Q"^ 

with q' = q - Pa, Q'^ = -q'^. 

(2,) 

To evaluate the contribution P^^ we rearrange the logarithmic term 

, Q'^ , (1 + a){Q' + XHa)) 
'"1? = • 

with 

^ ' 1 + a 

Performing the Borel transform and using the numerical values of the moments of the 
structure function t]^ (actually we used those, presented in [BS]) we find all the coefficients of 
the expansion to be of the same order of magnitude. Hence, following [30] we do not use 
this expansion here. We employ the explicit expression ^6], providing for the contribution to 
the LHS of Eq. (57) 

+ (4 + 5a)Gi{M^, a))r]a{a) - 18m^ar]b{a)]p. (76) 

Here G„(M^,a) = M'^'^"-'^^^ En{X'^ (a)) . We included only the lowest moments of the function 
rjb{o:), which are related to the moments of the structure function rja by expression, analogous 
to Eq. f l74|) - see |36]. The values of these moments enable to expect the convergence of the 
latter expansion. It will be useful (see next Subsection) to present the contribution in terms of 
those of lowest nonvanishing (second) moments of the structure functions. For the P structure 
it is 



ill 



3 m 

where 

"1 



A^2 = / daarjaia). (7J 
Jo 

The conventional notation is Ain = Iq daa'^~^'r]a{a). Numerically A^2 = 0.32 
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To avoid the cumbersome formula, we do not present the expressions for the contribution 
of the nonlocal vector condensate to the q structure of Eq. ljFri) . addressing the readers to our 
papers [3^ and [H]. The contribution of the second moment is 

(^4p)2 = -^mM^EoiM^, Wl)M2p. (79) 

Note that the nonlocality of the scalar condensate vanishes in the chiral limit [36], due to 
equations of motion. Thus B^^p = 0. 

The contribution of the condensate (M|gA^ ,j„G°'^'^cr^jyg|M) of dimension d = 5 vanishes in 
medium [26], as well as in vacuum [21] if we are using j = ji - see Eq. (|T3l) . There is no such 
cancellation for the currents j2 and J3. This is one more argument in favor of using the current 
ji, since the in- medium value of this condensate is not known. 



4.3 Approximate solution 

In this subsection we find an approximate solution of Eq. ( I57|) . The nucleon self-energies are 
expressed explicitly as functions of in-medium QCD condensates. 

This solution can be obtained by replacing the in-medium threshold value by its vacuum 
value in Eqs. (57). This can be done, since the variation of the threshold SW"^ contains a 
small factor e"^ . The accuracy of the approximate solution will be discussed in the next 
Subsection. 

We put = W"^ on the LHS of Eqs. f l59|) and multiply both numerators and denominators 
by g"^ 1^ _ Taking into account only the terms of dimension = 3 in the medium contributions 
we can write these equations as 

^^fp{M^)v{p) ^ ^ _ 

m*. (80) 



£g(M2)e™VM^ - g/,(M2)t;(p) 
£^(M2)e"^'/^^' - A'K^fi{M'')Kp{p) 



£g(M2)e-VM2 - 8^/^(M2)t;(p) 
Here the dependence of the contributions, provided by medium is contained in the functions 

,,,,2. [(3 - m')M^Eo{M') - M^E,{M^)]e^'l^' M'E,{M^)e^' l'^' 
Jq[^ ) = jjj-Q ; Jp[^ ) = jjj-g ; 

(81) 

fj{M^) = M^Ei(M2)e"^'/*^' 

One can see that they exhibit only weak dependence on in the interval (1271) - see Fig. 9. 
They can be approximated as 



/g(M2) ^ Co = 3.29; /p(M2) ^ cp = 1.30; fiiM^) ^ a = 1.60. 

Here the overline means that the function of is replaced by a constant value, corresponding 
to the lowest value of the function defined by Eq. (26) in the interval (27) at the vacuum 



9,9, 



value of the continuum threshold W"^. The RHS of these equations are written in powers of 
GeV (their dimensions are GeV^). 

Now one can replace the functions /j(M^) by the parameters q in Eqs. f lHOj) . Also, due to 
Eqs. (EOD and (EI]) one can replace the functions Cl{M^, W'^)e"''/^' and Ci{M\ W^)e"''/^^' by 
the constant values and mX^ correspondingly. We define also for the contributions induced 
by the matter 

. _ £^(M2,iy2))emVAf2 

Here the overline has the same meaning as in Eq. (I8ip . 
Now Eqs. ( I59l) can be written as 

Ey = : m* = —— . 82 

This is a direct analog of Eq. (50). We can identify 

= (83) 

We denote also 

= Y.J'n, (84) 

n 

where J-"^ is the contribution of the condensates with the dimension d = n. 

Employing Eq. (68) we find the contribution of the condensates with d = 3. In Eqs. (I82l) 
and (EH]) 

J-| = -A7vip); = 66i;(p); J"/ = -32«:,(p), (85) 

with all values written in powers of GeV. Using Eq. (3) we immediately find the values of 
and Sy at any value of the density p. For example, at phenomenological value of saturation 
density 

p = po = 0.16Fm-3 = 1.3-10"^GeV3 (86) 

we obtain 

Sy = 314 MeV, Eg = 0.18 . (87) 

Now we include the condensates of dimension d = 4. The contribution of the gluon con- 
densate can be obtained by employing Eq. f l69|) at = W^. We present the contributions of 
the nonlocal vector condensate at = through the contributions J-'f', corresponding to 
inclusion of only the second moments of the structure functions Ai2, given by Eqs. ( !77|) and 
(179|) . Direct calculation provides J-'^J^ = 290A^2P, and = 55AI2P in GeV units, while 

J-f = 0.37jf ^ ; J| = 0.62 . (88) 
Thus, with inclusion of the condensates of dimension d = 3,4 

7q = -47v{p) + 8.6gp + 34M2P • (89) 

Also, from Eq. (78) we find 

J-p = QQv{p) - 106A^2P , (90) 



while 

7i = -32Kp{p), (91) 

with the RHS written in GeV units. The value of J-"^ remained unchanged. Recall that A^2 = 
Jq daa7]a{a) = 0.32 is the second moment of the nucleon deep inelastic structure function. 



4.4 Gas approximation 

In the lowest order of OPE expansion the nucleon scalar self-energy depends on the scalar 
condensate K,{p) defined by Eq. (4). The next to leading order terms contain the gluon con- 
densate g{p), defined by Eq. fITOl) . Below we shall consider the problem of finding the density 
dependence of these condensates. In this Subsection we include them in the framework of the 
gas approximation [211 125] 

Kip) = K{0) + KMp; KM = {N\Y,q'q'\N), (92) 

i 

and 

9{p) = 9{0)+gNP; gN = {N\^G^'^^G';,jN). (93) 

71 

Thus in this approximation i^p = i^n Qp = Qn P- 

The scalar quark condensate is numerically more important. Numerous model calculations - 
see, e.g. [MIEH], and [70] for earlier papers, provide relatively small nonlinear term S{p) on the 
RHS of Eq. (4) at p close to the saturation value po- Thus, employing of the gas approximation 
is reasonable. 

Now we must find the values of k^t and g^- 



4.4.1 The value of {N\uu + dd\N) 

The expectation value is connected to the pion-nucleon sigma-term a by the relation [28] 

2(7 

kn = {N\uu + dd\N) = . (94) 

On the other hand, the a term can be expressed in terms of the pion-nucleon elastic scattering 
amplitude T = T{s,t,k'^,k''^) [28]. Here k and k' are the pion momenta before and after 
scattering, s and t are the Mandelstam variables s = (p + k)"^ and t = {k' — k)"^, with p denoting 
the momentum of the nucleon. The a term a is connected to the amplitude F in the unphysical 
point 

T(0, 0,0,0) = (95) 

J TT 

where /,r is the pion decay constant. The experiments provide the data on the physical ampli- 
tude Tphys = T{{m + ■mT,)^,ml,ml,m'^). The method of extrapolation of observable physical 
amplitude to the unphysical point was worked out in [22] • The S term H'^^ = —Tphygfl differs 
from the a term by 

= 15MeV. (96) 
Note that (S'^^ — o")/(T ~ m^r, in the chiral limit we must put S = cr. 



For many years the value 

a = (45±8)MeV, (97) 

based on the results of [71] was assumed to be true [72]. This corresponds to S'^^ fa 60MeV. 
However, some of the latest results support rather the value of S'^^ close to 80 MeV [TSj [TU [75] , 
corresponding rather to 

a = (60 - 70) MeV. (98) 

This is somewhat closer to some of results, claimed earlier [7B] . 

There is large discrepancy between the theoretical results - see [70] for references. Note 
that some of the latest ones [771 EH] support the value cr = 45 MeV. 

4.4.2 The value of {N\^G''f''G1;^jN) 

This expectation value was calculated in [79] by averaging the trace of the QCD energy- 
momentum tensor G over the nucleon state. The trace is 

= ^m,gg,-^G^ (99) 

where rrii the mass of the quark with the flavor i, G"^ = C^^^G^j,, 6=11 — 2n/3,nis the total 
number of flavors. There is a remarkable cancellation [79] 

{N\Y.m,m,\N) - ^1{N\^G'\N) = 0. (100) 

Here j stand for the "heavy" quarks with the masses rrij much larger that the inverse confine- 
ment radius /i (the accuracy of Eq. fllOOl) is {fi/rrijY)] is the number of heavy quarks. Hence, 
the charmed quark is the lightest one among those, which contribute to the LHS of Eq. fllOOp . 
and 

gN = -lim-Y^m,{N\Y^(tqW), (101) 
with i = u, d.s. The simplest expression for ^f^v can be obtained in the chiral S'f/(3) limit 

9N = g^N^ = -^-m. (102) 

One can see that 

Q>gN>g^N~^ . (103) 

Various estimations for the strange condensate in nucleon (A^|ss|A^) provide controversial results 
(see, e.g. [71]). However, since contribution of the gluon condensate is small even ioi gN = g^N i 
the estimation fll03p is sufficient for our analysis. 
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4.4.3 Numerical results in the gas approximation 

The results of previous subsection enable us to obtain the complete set of numerical results. 
We compare the approximate solution given by Eqs. f l85l) . f l89l) . fl90l) . fl9T]) and "exact" solution 
of the system of Eqs. (57) at the saturation value of density p = po- 

Start with the value kn = 8, which is consistent with Eq. fp7|) . Including only the conden- 
sates of dimension ci = 3 we find for approximate solution Sy = 314 MeV, m* —m = —205 MeV. 
Note that due to relatively large value Eg = 0.18 there is a contribution of about +200 MeV 
to m* — m, which is proportional to the vector condensate, and the one of about —400 MeV, 
proportional to the scalar condensate. Solving Eqs. (57) we find 

Sy = 323 MeV, m* - m = -201 MeV, 

- = -0.077 GeV^ Wl-W^ = 0.12 GeV^ (104) 

Thus, indeed the approximate solution is a very accurate one. 

Turning to the condensates of dimension d = 4 and using Eq. (11031) we see the gluon con- 
densate provides a small contribution of 0.009 to Eg in the approximate solution. Thus it adds 
about 3 MeV to both Sy and m*. The nonlocal vector condensate subtracts 0.014 from the 
value of Sq. The nonlocality of the vector condensate subtracts also 54 MeV from Sy. The 
approximate solution provides now Sy = 257 MeV, m* — m = —202 MeV. Inclusion of the 
condensates with dimension c? = 4 changes the solution of Eq. (157|) to 

Sy = 253MeV, m* - m = -213 MeV, 

A^^ - A^ = -0.23 GeV^ W^-W^ = 0.04 GeV^. (105) 

For K = 11, corresponding to a = 60 MeV fl98p we find 

Ey = 238 MeV, m* - m = -368 MeV. 

Fixing the continuum threshold = W"^ we would obtain 

Ey = 276 MeV, m* - m = -336 MeV. 

Thus the approximate solution, discussed above becomes less precise at larger values of k. 

We show the density dependence of the nucleon parameters and of the effective threshold 
value for kn = 8 in Fig. 10. The approximate solution provided by Eqs. fl89l) - fl9T|) is shown 
by dashed line in Fig. 10 a. In Fig. 11 we show the dependence of these parameters on at 
saturation value of nuclear density po- 

4.5 Asymmetric matter 

4.5.1 Condensates of the lowest dimension 

Now we employ the SR approach for calculation of the nucleon self-energies in the nuclear 
matter, composed of neutrons and protons, distributed with the different densities p„ and 
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Pp, i.e. in the SU{2) asymmetric matter. We calculate the dependence on the total density 
P = Pn + Pp and on asymmetry parameter 

/3 = (106) 

Pn + Pp 

We shall calculate the parameters for the proton. Those for the neutron can be obtained 
by changing /3 to — /3. As well as in the case of symmetric matter the lowest order OPE terms 
depend on the vector and scalar condensates. However, in addition to the isospin symmetric 
condensates v{p) and k{p), defined by Eqs. (3), (4), we have two SU (2) asymmetric condensates 
of dimension d = 3. These are the vector condensate 

vl^\p) = {M\u-f^u-d-f^d\M), (107) 

and a scalar one 

Cmip,l3) = {M\uu-dd\M). (108) 

The lower index m is introduced in order to show that for the SU{2) invariant vacuum Cm(0) =0. 

The vector isotope asymmetric condensate can be immediately calculated. In the rest frame 
of the matter 

ti-\p) = vir\p,P)K; v^-\p.P) = -Pp. (109) 

The situation with the scalar condensate is more complicated. Even in the gas approxima- 
tion, where 

Cm(p, /3) = -/3pCp; Cp = {p\uu - dd\p) (110) 

cannot be directly related to observables. Thus we need some model assumptions. 

Considering the nucleon as a system of three valence quarks the sea of the quark-antiquark 
pairs we can write 

Cp ~ Cp + Cp ) 

where the upper indices denote the contributions of the valence and sea quarks correspond- 
ingly. In nonrelativistic quark models = 1, while Cp < 1 for the relativistic models. The 
experimental data on the deep inelastic scattering provide Cp ~ —0.15 [80l [8T] . 

4.5.2 Nucleon self-energies 

In this Subsection we include only the condensates of the lowest dimension d = 3. We can 
write the contributions to the LHS in a form similar to Eqs. ( 1671) . The term A\p depends only 
on the vector condensate f (p), and does not depend on j3. Thus it remains the same as for the 
symmetric matter, while 

= -47rVEi(M2, ^2)^^^^^). 



Ap=-^ ^^"^v^:;"^-^ ^p,/3)- (in: 



ill- AM^Ei{M^,W^ 
'3 

The functions £"„ are given by Eq. (23). Here we defined 9{p, (3) = (k^v + PCn)p and w(p, (3) = 
3p(l — /3/4). Comparing these expressions with Eq.( l671) we find a simple solution in which the 
proton self energies are expressed in terms of the solution for symmetric matter [33] . 

Ey(p, /3) = Sy(p, 0)(1 - /3/4); m*(p, k^, Cp) = m*{p, + /3Cp, 0, 0); (112) 
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Wl{p,P) = Wl{p,Q). 

For example, employing Eq.f ll04p we find that in the neutron matter (/3 = 1) the proton 
vector self-energy is 242 MeV, while for neutron it is 403 MeV. We shall see that inclusion of 
the higher condensates will modify the values. 



4.6 Possible mechanism of saturation 
4.6.1 A special role of the pion cloud 

Note that the simplest account of the nonlinear effects signals on a possible saturation mech- 
anism in our approach [2l]-[26]. The nucleons of the matter interact by meson exchange, and 
the nonlinear contribution to the scalar condensate S{p) - see Eq. (4) comes from the meson 
cloud, created by nucleons. The contributions of each meson is thus proportional to the value 
of the operator qq averaged over the meson states. 

We expect exchanges by pions to dominate the value S{p). The matrix element {h\ 
counts the total number riqq of quarks and antiquarks in a hadron h ESI IM] under certain 
reasonable assumptions on the quark wave function of the hadron [831 El] • Thus one can expect 
{pj\(pq^\fij) ~ 2 for a meson fi. However for the pion, which is a collective Goldstone excitation 

2m^ 

{7c\Y.Q'Q» = ^ = 2m^^5, , (113) 

(the factor 2m^ comes from normalization of the pion vector of state). Thus riqq ^ 12 for pions, 
and we expect their exchanges to provide the leading contribution to S{p). 



4.6.2 Lowest order in Fermi momentum series 

Here we consider the case of symmetric matter. We shall calculate the contribution the non- 
linear term S{p) in the lowest order of expansion in powers of Fermi momentum. The latter is 
related to the density as 

P = ^P%- (114) 

The leading contribution is provided by single-pion exchange term (known also as the Fock 
term or the Pauli blocking term). It can be written as 

S{p) = -j ^^0^mD^'^\k)TmAuu + dd\7^). (115) 

where pi^2 ^ Pf are the nucleon momenta, k = Pi — P2, T stand for the irNN vertices, 
= 1/(0;^ — k'^ — + ie) is the propagator of free pion carrying the energy u. Summation 
over spin and isospin variables is assumed. The corresponding contribution to the correlator is 
illustrated by Fig. 12. The function S{p) is obtained analytically [36]. The expression becomes 
especially simple in the chiral limit = 

Q( \ 9 PFp2ml 
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where the last factor is the pion expectation value Eq. flll3p . Thus S{p) ~ p^^^. Note that 
the heavier mesons do not contribute to the leading term of the Fermi momentum expansion. 
Their contribution is of the order p^. 

The difference between the RHS of Eq. flll6p and the result, obtained with the use of the 
physical value of is rather large. However, in the chiral limit one should use the value of S 
term, rather than that of cr-term for the calculation of k^- - see Eq. (96), since the difference 
E — 0" contains additional powers of m^. This diminishes the difference. As a result, the value 
of K with kat replaced by 

< = (117) 

appears to be close to that, obtained for the physical pion mass [30]. Some additional arguments 
supporting the use of the chiral limit at p close to po were given in [85]. 

Now we expand the nucleon potential energy f/(p), defined by Eq. (55) in powers of Fermi 
momentum. We take into account the contributions of the order p and ppi? ~ p^^^, neglecting, 
however, the terms of the order p^. Including only the condensates of the lowest dimension 
d = 3 and using Eq. (|85|) we find (in GeV units) 

= 66v{p); = A6v{p) - 32Kp{p). (118) 

Here Kp{p) = K^p + S{p) with S{p) given by Eq. (11161) . Thus the nucleon energy on the Fermi 
surface can be written as 

p{z) = [38^2/^ + (430 - 42Kfl)z + 126^/^] MeV. (119) 

Here z = p/po, with po being the phenomenological value of saturation density - Eq. ( 186|) . The 
first term on the RHS is the kinetic energy. 

4.6.3 Equation of state 

The saturation value p can be found by minimization of the binding energy per particle 

e(p) = - r dxp{x) . (120) 
p Jo 

In our case 

e{z) = [23z^/^ + (215 - 21k^)z + 54^^/^] MeV. (121) 

One can see that e(z)' = at z = 1 if nf^ = 14.4. This corresponds to S'^^ = 79 MeV, in 
agreement with the experimental data. The nucleon energy at the Fermi surface is p = —11 
MeV. The vector self-energy is ^ 260 MeV, while m — m* ^ —310 MeV. The value of 
incompressibility K = Qp^d^e/ dp^ = 170 MeV. A standard value obtained in the nuclear physics 
approaches is K ^ 230 MeV [86]. Inclusion of the condensates with d = A modifies these values 
slightly. 

Note that somewhat similar saturation mechanism take place in the model, suggested later in 
[87] . The approach of [87] was based on the chiral effective Lagrangian. The authors calculated 
the nonlinear contribution to the scalar condensate, caused by the pion exchange. The value 
of S{p) appeared to be close to the ours one. 
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Of course, one should not take these results too seriously. We tried to obtain the values of 
the self-energies. The calculations of the potential energy require somewhat higher accuracy, 
since it is a result of subtraction of two large values. Thus, reasonable values for the potential 
energy is a surprise. Of course, we can not expect accurate values for the binding energy. 
Also, the results are very sensitive to the exact value of k^. For example, in this approach 
the saturation is achieved at p = 2po if i^i} = 14.7, corresponding to S = 83MeV. Also, the 
contribution of two-pion exchange to S{p) is too large to be neglected [5S] . 

However, the results of this subsection can be treated as a sigh of a possible mechanism for 
saturation of nuclear matter. It is due to the nonlinear behavior of the scalar condensate. Thus 
in our approach the saturation is possible only if the three-body interactions (strictly speaking, 
many-body interactions) are included. 

The potentials of the form U{z) = Ciz + 6*4/32;''/^ were studied earlier in nuclear physics. 
Such potential with Ci ^ — 210MeV, C3/4 160 MeV was analyzed in [89]. Note that our 
values in Eq.([Il9]) are Ci = -175 MeV (for = 14.4 )and C4/3 = 126 MeV. 

Note that this mechanism differs from that of the Walecka model [H [2] . Recall that in the 
latter case the saturation is due to the different behavior of the vector and scalar fields on 
density. The vector and scalar fields can be expressed as [I] 



with the sum over spin and isospin variables being assumed. Here and gs are the coupling 
constants. In the vector case = un{p)joUn{p)/'2'E{p) with un standing for the nucleon 
bispinor, while E(p) = V{p) + {m*(p) + p^)^/^. One finds immediately that A^^, = 1, and 
V{p) = g^p is exactly proportional to p. In the scalar case Ns{p) = uj\f{p)uj\f{p)/2E{p), leading 
to a more complicated dependence $(ps) = QsPsip): where the scalar density 



is a nonlinear function of p. Thus, the saturation in Walecka model is a relativistic effect. It 
takes place in framework of the two-body interactions. On the contrary, in our approach it does 
not require relativistic treatment of the nucleons of the matter, being caused by the many-body 
interactions. 

We shall return to this problem in Sec. 9. 

5 Other characteristics of the in-medium nucleons 
5.1 Axial coupling constant 

The SR method have been used for calculation of the nucleon coupling constants for isovector 
and isoscalar axial currents in vacuum |90]-[92]. The approach is based on considering the 
function Hq defined by Eq. (fT2l) and the dispersion relation ([9]) in external axial field A^, coupled 
to the quarks by interaction V = gqA^j'^'j^, and included in the lowest order of the perturbation 
theory. In the isovector case, directly related to the neutron /3 decay = l,gd = —1. The 
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corresponding nucleon coupling constant in vacuum is qa = 1-27 Recall that the vector 
coupling constant is = 1, and the nonzero difference qa — Qv is due to the nonconservation 
of the chiral current. 

The function Ho contains several structures, which are linear in A. The structure q{A- 
g)75. appeared to be the most convenient one for the SR analysis. The LHS contains now 
several expectation values, which vanish in the absence of the axial field. The condensate 
'•a = {^\<iG ai3l i3(l\^) A = ^ai^(0 IggjO) , with V usually referred to as the field induced susceptibility 
appeared to be numerically the most important one. The numerical value of this expectation 
value was obtained in (93j . 

On the RHS of the SR the lowest lying singularity is now the double pole, corresponding to 
the contribution (0|j|A^)(A^|A75|A^)(A^|j|0). The SR reproduced the experimental value of the 
isovector coupling constant qa ~ 1.25. 

The first attempt to calculate the value of q'a for the nucleon in nucleon matter was made in 
j94j . The authors employed Eq. (30), combining it with several phenomenological assumptions. 

In [38] the renormalization of the nucleon coupling constant Sqa = 9^ — 9 a in nuclear matter 
was calculated by extending of the method developed in [23]-[25] for the case of the external 
axial field. The LHS of the SR was dominated by the configuration, in which one of the qq pairs 
was exchanged with the condensate ta-, and the other one-with the quarks of the nucleons of 
the matter. The RHS of the SR included the contribution of the isobar-hole (AN) excitations 
in terms of the single-pole terms. 

The result of calculation depend on the value of the scalar condensate in nuclear matter. 
Thus it depends on the actual value of the a term. It was obtained that 6gA = —0.24 at the 
saturation density if cr = 60 MeV. For a = 45 MeV the shift is Sqa = —0.22. This is in agrement 
with the experimental data [951 ES], providing g^{po) ~ 1.0. If the single-pole terms are not 
included, the result is Sqa = —0.05. Hence, the A — iV excitations is the main mechanism of 
the process, in agreement with [9| [TO]. 

However, the approximations for the LHS of the SR, made in [3S] are too crude. A more 
rigorous analysis, which includes the four-quark condensates is still needed. 

5.2 Charge- symmetry breaking forces 
5.2.1 Neutron— proton mass difference 

It is known [56] that the difference between the neutron and proton masses in vacuum is m„p = 

1.3 MeV, while the electromagnetic interactions contribute to this value is m^^ ~ —0.7 MeV [H] . 
Thus the contributions of the strong interactions is m^^ ~ 2.0 MeV. If the isospin symmetry 
(known also as the charge symmetry) in the hadron interactions is assumed, the value of m^^ 
would not change in nuclear matter. Thus, it is interesting to study the density dependence 

Kp(p)- 

One can expect the mass difference of mirror nuclei to be 

AM = E'^ + mnp, (124) 

where is the electromagnetic energy difference. The values of AM and E'^ can be, corre- 
spondingly, measured and calculated with high accuracy [97]. However, this equation is not 



satisfied if m„p is equaled to its vacuum value [OB]- The discrepancy between the two sides of 
Eq.(124) increases with the atomic number A, reaching 0.9 KeV for A = 208. This is called 
the Nolen-Schiffer anomaly (NSA). 

The phenomenological charge-symmetry breaking (CSB) potentials of the nucleon interac- 
tions were discussed in [11]. Some of them described the NSA, but contradicted the data on the 
CSB effects in NN scattering. A possible resolution of the NSA is the assumption that the value 
of rrinp in a nuclei is smaller than that in vacuum. In the case of nuclear matter this assumption 
would manifest itself in the decrease of mnp{p) while p increases. The renormalization of m^p 
was calculated in numerous papers (see, e.g. [99]- |103j ) for nuclear matter and for the finite 
nuclei. However, the NSA is not fully understood yet. 

5.2.2 QCD sum rules view 

The QCD SR provide a consistent formalism for the CSB effects. In QCD the charge-symmetry 
is broken due to a nonzero value of the difference of current quark masses 

rud ^ 7MeV; m„ ^ 4MeV; /i = - m„ ^ 3MeV. (125) 

Besides the quark mass difference the CSB manifests itself also through a nonzero expectation 
value of the operator uu — dd. In vacuum the ratio 

{0\dd-uu\0) 
^° ~ {0\uu\0) 

was obtained by the QCD SR approach [104[ I105[ 1106] basing on the experimental values of 
isospin breaking mass splitting of nucleons and hyperons. 

In nuclear matter one can expect the charge-symmetry breaking effects to be expressed in 
terms of the quark mass difference /i and of the isospin symmetry breaking condensate 

{M\dd-uu\M) 
^"^ ~ {M\uu\M) 

The first attempt to calculate mnp{p) was carried out in |107] . The analysis was based on 
the vacuum SR, obtained in |104] . Only the terms, containing scalar condensate were included. 
The vacuum scalar condensate in equations of |104j was replaced by its in-medium value, 
obtained under certain phenomenological assumptions. Calculations of |108[ 1109] included also 
the density dependence of the CSB parameter 7^. 

The SR calculations of [39] present the direct extension of the approach, developed in [36] . 
Now we must calculate the contribution of the strong interactions to the difference of the 
binding energies e„p = En — Sp of the neutron and proton in symmetric nuclear matter. One 
can write Ei = Ui + Ti, with Ui{Ti)- the potential (kinetic) energy of the nucleon i = p,n. The 
potential energies are expressed in terms of the self-energies by Eq. (55), while T, = p|,/2m*. 

Now the LHS of Eqs. (57) should be calculated for the neutron and proton separately. The 
explicit dependence of the free quark propagators gq, defined by Eq. (15) on the quark masses 
should be taken into account. Also, the term £^ on the LHS of Eq. (57) contains the CSB 
condensate {M\uu — dd\M). Thus the difference of the binding energies can be written as 

Enp{p) = phi{p) + '^rn{p)h2{p) ■ (126) 



The functions &i,2(p) can be obtained by the SR method. The density dependence of the 
condensate 7^ can be obtained in framework of certain models. 

Start with the calculation of 61,2 The terms, proportional to the quark mass difference 
contribute to the difference of the scalar self-energies 

5, =0.17^^, 
Po 

and to the difference of the vector ones 

^2 = -0.029^^. (128) 
Po 

Recall that v{p) = 3p is the vector condensate, Hp{p) = n{p) — /t(0), with the scalar condensate 
K defined by Eq. (4), po is the value of saturation density, given by Eq.( !86l) . There is also a 
contribution to the difference of the scalar self-energies, containing the CSB condensate 7^ 

^3 = C^(7n^('^(p) - + (7™ - 7o)«:(0)); C = 32 GeV-2. (129) 

One should include also the strong-interaction part to the neutron-proton difference of the 
vacuum parameters m„p, A^^ and W^^. The two last ones were obtained in |105] by the SR 
method, basing on the empirical value of m'^np- These parameters, as well as the CSB condensate 
7o were expressed through the quark mass difference p. 

We present the values of 61^2 at the saturation value of the density 

&i(po) = -0.73; 62(po) = -l.OGeV. (130) 

The density behavior 7m(p) = 7o('^(p)/'^(0))^^^, based on the Nambu-Jona-Lasinio model 
was suggested in |108] . Under this assumption and employing 70 = — 2 ■ 10"'^ |105] we obtain 
£^np(po) = — 0.9MeV. Our final result is 

^np(po) = (-0.9 ± 0.6) MeV, 

with the errors caused mostly by uncertainties of the value of 70. 
Thus at least qualitative explanation of the NSA is achieved. 

5.2.3 Consequence for conventional nuclear physics 

At least two points of the QCD SR analysis can be useful for nuclear physics. 

In conventional nuclear physics the u — p mixing in the vector channel is usually believed 
to be responsible for the largest part of the Nolen-Schiffer anomaly . However, the SR analysis 
predicts the scalar channel to be important as well. Neglecting the terms, containing the scalar 
condensate, we would obtain Enpipo) > 0, This contradicts the experimental data and general 
theoretical expectations. 

The other point is the Lorentz structure of the nucleon interactions. In the QHD Eq. ([1]) 
the terms V and $ are caused by the vector and scalar interactions correspondingly. In the SR 
analysis described in Sec. 4. these terms are determined by the vector and scalar condensates 
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correspondingly. Such separation is violated by inclusion of the mass terms in the quark prop- 
agators Qq determined by Eq. ( IT^ . Relation between the Lorentz structures of Eq. ([T]) and the 
condensates becomes more complicated. Now there are contributions to $, which are propor- 
tional to the vector condensate and contributions to V , proportional to the scalar condensate. 
These admixtures are small, being proportional to the current quark masses. This anomalous 
Lorenz structure are shown in Fig. 13. 

In the case of the CSB forces they manifest themselves in the leading terms. One can see 
that 5i given by Eq. (11271) contributes to the structure $ of Eq. (1), being determined by the 
vector interaction. On the other hand, the contribution 82 presented by Eq. (I128P is determined 
by interaction in the scalar channel, but contributes to the vector structure V . This correspond 
to anomalous structure of the nucleon vertex functions. 

These points can be helpful in constructing the CSB nucleon forces. 

5.3 Nucleon deep inelastic structure functions 

The SR approach enables to investigate internal structure of the nucleon. Such problems are 
unaccessible for nuclear physics. 

5.3.1 Nucleon structure functions in QCD sum rules 

In deep inelastic scattering (DIS) the electrons transfer large energies and momenta to the 
target. This process is a well known tool for investigation of the internal structure of the latter. 
The cross section of DIS can be expressed through the imaginary part of the amplitude of the 
elastic scattering of the virtual photon. The latter has a large and negative four-momentum 
squared /c^ < 0, —k"^ ^ m? (m is the rest mass of the target) and the energy ^ m, 
while —k'^/ko ~ m. The investigation of DIS enables to study the momentum distribution of 
the quarks. The structure functions F2{x) measured in the DIS are determined by the quark 
distributions xq{x), where x is the nucleon momentum fraction carried by a quark. 

The QCD SR method was applied for investigation of the DIS on the proton. The second 
moments of the structure functions were obtained in |11U] and |111] . The next to leading terms 
of the asymptotic expansion in powers of were found in [112j . The structure functions 
F2{x) at moderate values of x were calculated in pT3]. Another presentation of the structure 
functions was obtained in |114j . We shall employ the approach, developed in jll4j since it can be 
extended for the case of finite density in a natural way. On the other hand, such generalization 
is the extension of approach described above. 

In order to obtain the distribution of the valence quarks, the authors of [114] considered the 
correlation function G , which describes the system with the quantum numbers of the proton, 
interacting twice with strongly virtual hard photons 



Here q and q + k are the momenta of the nucleon in the initial and final states, k = ki — k2 is 
the momentum, transferred to the system by the photon scattering. The incoming (outgoing) 
photon carries momentum ki{k2), interacting with the quarks in the point y — A/2 {y + A/2). 
The quark-photon interaction is described by the function H{y,A). 
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In |114] the correlation function G{q^ k) was calculated in terms of the QCD condensates. 
The double dispersion relation in variables = and = {q + k)"^ was employed. The 
crucial point was the OPE in terms of nonlocal operators depending on the light-like vector A 
(A^ = 0). The Borel transform in qf and gf was carried out. The equal Borel masses = M| 
are considered. The Fourier transform in A provided the momentum distribution of the valence 
quarks. 

Thus, in |114] the momentum distributions of the valence quarks were expressed in terms 
of the QCD condensates. 

Before discussing the extension for nuclear matter, we describe the most intriguing problem 
of the DIS on nucleus. 

5.3.2 EMC effect 

The experimental data obtained first by the EMC collaboration |115] showed that the DIS 
functioning"^ (x) of nucleus with the atomic number A differs from the sum of those of free 
nucleons. The structure function was compared to that of deutron [A = 2), which imitates the 
system of free nucleons. The deviation of the ratio 



from unity characterizes the deviation of a nucleus from the system of free nucleons. 

The ratio R"^{x) appeared to be the function of x indeed. Most of the early data were 
obtained for iron (Fe). Exceeding unity at x < 0.2, the ratio drops at large x, reaching 
the minimum value R^^ ~ 0.85 at a; ~ 0.7. This behavior was called EMC effect. The 
same tendency was traced for other nuclei. Both experimental [116] and theoretical jll7] 
investigations of the effect are going on nowadays. 

There are several mechanisms which may cause the deviation of R^{x) from unity. We shall 
try to find how the difference of the quark distributions inside the in-medium and free nucleon 
changes the ratio R^{x). 

5.3.3 EMC effect in QCD sum rules 

In [lO] the approach of [114j was combined with that of [36] for calculation of the quark 
distributions in the proton, placed into the nuclear matter. The expression for the correlation 
function took the form of Eq. (11321) . with the averaging over vacuum replaced by that over the 
ground state of the matter. The two types of contributions to the correlator were considered 
- Fig. 14a,b. In the diagram of Fig. 14a. the photon interacted with the quark of the free 
loop. In the diagram of Fig. 14b. it interacts with the quark, exchanging with the matter. The 
modification of the distribution of the quarks was expressed in terms of the vector condensate, 
with its nonlocal structure being included, and through the shift of the scalar condensate. The 
results are true only for moderate values of x. They can not be extended to the region x ^ 1, 
since the OPE diverges in that region |113] . 

Omitting the details of calculation, provided in [lOj, we present the results in Fig. 15. We 
carry out calculations for nuclear matter. The ratio R{x) can be viewed as the limiting value 
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of R^{x) at A — i- oo. One can see that the distributions of u and d quarks in fraction of the 
target momentum x are modified in a different way. The fraction of the momentum carried by 
u quarks decreases by about 4%. The ratio R, determined by Eq. fll33p has a typical EMC 
shape. 

Note that while for a free nucleon x < 1, for the nucleus it can be as large a.s x = A. The 
developed approach enables to calculate the quark distributions at x > 1, describing thus the 
cumulative aspects of the problem. 

5.3.4 Swelling of the nucleon 

One of the possible explanations of the EMC effect is the "swelling" of the nucleon inside the 
nucleus. Under this assumption the structure function of the in-medium nucleon F^{x) is 
written in terms of that for a free nucleon F2{x) as F^{x) = (xr^/To); with < tq |118] . 
On the other hand, the value of the proton residue is proportional to the square of the 
three-quark wave function at the origin. Assuming that there is only one length scale r for the 
nucleon, we find ~ r~^. Employing Eq. (11051) we find (r^ — 'ro)/rQ = 2%, in agreement with 
estimations of |119] . 

6 Intermediate summary 
6.1 Reasons for optimism 

One can see that we have some reasons for optimism. The QCD sum rules analysis confirmed 
that the nucleon in nuclear matter can be treated as moving in superposition of strong vector 
and scalar fields of the order of about 300 MeV. The vector field is positive, while the scalar field 
is negative, and a partial cancelation takes place. We obtained this result without employing a 
controversial conception of interaction of the point-like nucleons. 

We obtained a picture for formation of the nucleon self energies. The effective mass m* is 
formed by the exchange of noninter acting quarks between the system of the three noninteracting 
quarks described by the current ji determined by Eq. (|T3|) which carries the proton quantum 
numbers and the scalar condensate, which differs from that in vacuum. Of course, there is 
strong interactions between the quarks, which form the condensate. The vector self energy is 
formed as the exchange by noninteracting quarks between the system of the three noninteracting 
quarks, carrying the proton quantum numbers and the valence quarks of the matter, forming the 
vector condensate. The quarks, which form the condensate interact strongly between themselves 
and with the other quarks of the matter. 

The vector and scalar self-energies are calculated in terms of the in-medium vector and 
scalar condensates, the nonlocality of the vector condensate is included. The vector condensate 
can be calculated easily. The contributions, corresponding to the nonlocal structure of the 
vector condensate can be presented in terms of the moments of the structure functions. Thus, 
it is also related to observables. The in-medium scalar condensate in the gas approximation, 
which has a good accuracy near the saturation point can be related to the observable pion- 
nucleon S term. Hence, the nucleon self-energies are expressed in terms of the condensates, 
which can be either calculated in a model-independent way or related to observables. 



The approach was used for calculation of other nucleon parameters, providing reasonable 
results. 

6.2 Reasons for scepticism 

However, the results, described above leave some reasons for scepticism. The most important 
one is the obscure role of the four-quark (4g) condensates. A simple estimation for the value of 
the scalar 4q condensate is {M\qqqq\M) = 2(0|gg|0)(A^|gg|A^)p. Employing this estimation one 
would find both vector and scalar self-energies to be large (hundreds MeV) and positive. This 
would contradict to the QHD phenomenology. 

Another obscure point is the role of the radiative corrections. As we have seen in Sec. 2, in 
the lowest order radiative correction to the vacuum SR the coupling constant is multiplied 
by a numerically large coefficient. The role of radiative corrections in the vacuum case was 
clarified in [M]. However, the case of finite densities require a separate analysis. 

Some of the results obtained in |120j are sometimes viewed as the grounds for additional 
sceptical statements. One of them is a possible strong shift of the hadron mass rrih due to 
the resonances in the Nh system, which is not included to the SR analysis. Another one is 
about the importance of the large distances for the formation of the nucleon mass, while the 
SR actually deal with small distances of the order 1 GeV~^. 

6.3 Response to the sceptical remarks 

Start with the latter statement. The nucleon wave function is indeed formed at large distances 
from the center of nuclei. The case of deuteron is a bright example. However, the NN inter- 
actions take place at much smaller distances. One can obtain a rather good description of the 
deuteron by employing the potential V{r) = C6{r). In Walecka model [1] the NN interaction 
radius is of the order of co meson inverse mass l/rrii^. The duality interval (27), where the SR 
are used, just corresponds to the distances responsible for the NN interactions. Hence, the SR 
approach should provide an adequate description of the nucleon in nuclear matter. 

Resonances in Nh system are described by singularities in variable s of the function Umiq"^, s), 
defined by Eq. (!38|) . We avoided these singularities by considering Ilm{q^, s) at fixed value of s. 

Contribution of the 4g condensates and those of the radiative corrections should be included 
into our analysis. This will be done in next Sections. Note that the calculations of the scalar 
4g condensate in nucleon in framework of the Nambu-Jona-Lasinio model [121J provided en- 
couraging results. The estimation, which we gave in the beginning of this subsection overshoot 
the value of condensate since there are certain cancelations. 



7 Four-quark condensates 



7.1 General equations for contribution of 
the four-quark condensates 

We shall need the expectation values 

T^yjuh = (M| : g^i^r^g^i'^'g^^^r^g^^^' : \M){6aa'Sw - ^a^'^^a'), (134) 

where the colon signs denote the normal ordering of quark operators, /1.2 stand for the quark 
flavors, a, a', b, b' are the color indices. A nonzero contribution to the function Um is provided 
by the antisymmetric combination of colors - see Eq. (1131) . The basic 4x4 matrices F"^'^, 
acting on the Lorentz indices of the quark operators are 

r^ = /; r^^ = 75; 1^ = 7^; ^^ = 1,15; = '-{i,iu-i,iu). (iss) 

Thus, they describe the scalar, pseudoscalar, vector, pseudovector (axial) and tensor cases 
correspondingly. 

In order to simplify the formulas, we shall not display the color indices, keeping in mind 
that the quark operators are color-antisymmetric. 

The contribution of the four-quark condensates to the LHS of the sum rules is |32] 

(n^ - no)4, = (n,)4, = 4 ( E f^xvH'''' + ^ r^y i?^"") . (136) 

1 X,Y X,Y 

Here 

H""^ = {M\uT''uuT^u\M); i?^^ = {M\dT'' duT^u\M), (137) 

while 

I^XY = ^Tr(7«r^7^r^)757"g7^5 ; 

TxY = ^Tr(7,g7^r^)757"r^7''75 • (138) 

The factor = 1 if has a vector or tensor structure, in other cases 6y = —1. The sign is 
determined by that of the commutator between the charge conjugation matrix and Fy. 

Considering the iu condensates, one can see that the contributions [IxyH^^ to 11^'^ obtain 
nonzero values only if the matrices Tx and Fy have the same Lorentz structure. All the struc- 
tures, presented by Eq. (I135P contribute to the Au condensate. In the case of 2d2u condensate 
the factor Tr(7Q,g7^F^) does not turn to zero only if Fy has a vector or axial structure. In the 
latter case Fx should have axial structure as well. In the former case it can be either Lorentz 
scalar or Lorentz vector. 

7.2 Approximations for the four-quark condensates 

We shall make certain approximations for the quark condensates both in vacuum and in nuclear 
medium. 



7.2.1 Factorization of the vacuum condensate 

The vacuum 4q condensate will be considered in framework of the factorization hypothesis [20] 

{0\uT^uuT^u\0) = {0\uT^u\0){0\uT^u\0). (139) 

Thus, on this approximation the condensate {0\ur^ uuT'^ u\0) does not vanish only in the scalar 
case = = /. 

{0\uuuu\0) = {0\uu\0){0\uu\0). (140) 

This approximation has been justified in the limit of large number of colors [22] • This is a 
standard approximation in the SR calculations. As it stands now, there are no indications of 
noticeable violation of the factorization relations Eqs. f ll39p . f ll40p . 

7.2.2 Gas approximation 

The modification of the 4g condensates in nuclear matter will be treated in the gas approxima- 
tion. We shall go beyond in Subsection 7.5. In the gas approximation expectation value of any 
operator X is 

{M\X\M) = {0\X\0) + p{N\X\N) , (141) 

where 

{N\X\N) = J d^x[{N\X{x)\N) - {0\X{x)\0)), (142) 

is the excess of the density X{x) over the vacuum value inside the nucleon. Of course, (0|X(a;) |0) 
does not depend on x. 

Assuming that at any space point x 

{N\uT^uuT^u\N) = {{N\u{x)T^u{x)\N)y, 

one can write Eq. 01421) for 4u condensate in another way 

{N\uT^uuT^u\N) = J d^x[{N\u{x)T^u{x)\N) - 

- {0\u{x)T^u{x)\0)y + 2{0\uT^u\0){N\uT^u\N) + 

+ 1/jv(((0|mm|0))^ - {0\uuuu\0)). (143) 

Here Vn is the volume of the nucleon. The last term vanishes under the vacuum factorization 
assumption - see Eq. (11401) . For all structures but the scalar one the second term vanishes, 
since {0\uT^u\0) = 0, and Eq. fll43p takes the form 

{N\uT^uuT^u\N) = J d^x[{N\u{x)T^u{x)\N) - (0|u(a;)^^M(x)|0))^ (144) 

However, for the scalar case Tx = / we can write 

{N\uuuu\N) = J d^x(^{N\u{x)u{x)\N) - (0|mm|0))^ + 2{Q\uu\0){N\uu\N). (145) 



In similar way we can write for the 2d2u scalar- vector condensate, which contributes to the 
RHS of Eq. (^SB forrcjy 



{N\ddu^ou\N) = J d^x[{N\d{x)d{x)\N) - {0\dd\0)) x 

X {N\u{x)-fou{x)\N) + 2{0\dd\0) {N\u-fou\N) . (146) 

In the gas approximation we can write for the condensates, defined by Eq. fll37p 

H""^ = h^'^pp + h^'^pn, h^'^ = {N\uT''uuT''u\N) (147) 

R^y = r'^'^p- r^^ = {N\dT''duT''u\N). 
Further calculations require certain quark model of nucleon. 

7.3 Perturbative Chiral Quark Model 
7.3.1 Description of the model 

The Perturbative Chiral Quark Model (PCQM) originally suggested in [30], was developed 
later in [122j . Further applications are reviewed in |123j . The nucleon is considered as a 
system of three relativistic valence quarks, moving in an effective static field. The valence 
quarks are supplemented by a cloud of pseudoscalar mesons, introduced in agreement with 
the requirements of the chiral symmetry. In the SU{2) version of the model, which will be 
used here, only the pions are included. The meson cloud is included in the lowest order of 
perturbation theory. 

The PCQM Lagrangian is the sum of the terms, describing the constituent quarks, pions 
and their interaction 

L = Lq + L^ + Unu Lq = i,{r)[id - S{r) - 7oF(r)]^(r); 

= \{dMr)f: Unt = -^V'(r)5(r)75^^^^(r). (148) 

Here 0j is the isotriplet of the pion fields. The contribution L^nt is the lowest order expansion 
of the chiral interaction term Lint = —'ip{r)S{r)e'^""=''^'''^^^^^^^^ip{r), /^r is the pion decay constant. 
We did not write down the quark mass term. 

Important steps in development of the model were made in |124j . In earlier applications the 
pions were considered as independent point-like degrees of freedom. In |124] they were treated 
as containing the "sea" quarks [qq pairs) of the nucleon. 

Another move touched the treatment of the constituent quarks. The authors of [124J did 
not solve the Dirac equation for given form of the potentials S{r) and V{r), but postulated the 
Gaussian shape of the constituent quark density. It was assumed that the coordinate part of 
wave function of the constituent quark can be represented as the product of three single-particle 
functions 

^(r) = <P{r)x{r); <P{r) = Ne-^'/'^'; x{r) = ( ^^k) ] (149) 

4n 



with the normahzation condition / d?ril){x)'-)Qip{x) = 1 expressing the conservation of the baryon 
charge. Parameters of the model /3 and R are fitted to reproduce the values of the axial coupling 
constant and of the proton charge radius correspondingly. Employing the Dirac equation one 
finds that the wave function represented by Eq. (11491) corresponds to the scalar field 

S{r) = M + cr^, 

where M = (1 — 3/3^)/2/3i? ^ 230 MeV can be treated as constituent mass of the quark. 
7.4 Four-quark condensates in the PCQM 

Now each of the condensates h^"^ and r"^^ defined by Eq. (I147P can be presented as consisting of 
three contributions. All the four operators can act on the valence quarks. This will be denoted 
by the lower index val. All the four operators can act on the pions. This will be denoted by 
the lower index P. There is also a possibility that two operators act on the constituent quarks 
while the other two act on pions. Such interference terms will be denoted by the lower index 
J. Hence, we can write 

h§ = {h§U + {h'')p + {h'')j- r^^ = {r^'^Ui + (r^^)p + (r^'')j. (150) 

Since in the case of 4-u condensate only the structures with X = Y contribute to the SR, we 
used notation h^-^ = . In the analysis presented in the next subsubsection we omit several 
terms, which are numerically not important. Start with the contribution of the valence quarks. 

7.4.1 Contribution of the valence quarks 

Here the products of the quark operators are averaged over the state |A^) of three constituent 
quarks, described by the wave functions presented by Eq. (I149p . Due to the normal ordering 
of the quark operators we find for neutrons {h^)vai = 0, while for protons 

{h^Ui = {U\uT''uuT''u\U). 

Here \U) denotes the vector of state of two constituent U quarks. 

For the scalar case the product of the PCQM operators describes the excess of the quark 
density with respect to the vacuum value 

uPCQM ^r^^PCQM ^ M(r)n(r) - {0\uu\0); 

(^jj^^PCQM^PCQM^^^ ^2 J d^rij{r)^{r). (151) 

Thus, employing Eqs. (11421) . (11431) we find for the scalar case X = Y = S 

{h^^ai = 2{0\uu\0){N\uu\N) + J d^r{i>{r)ij{r)y. (152) 

The last factor (NluulN) in the first term on the RHS has the meaning of the contribution 
of the valence quarks to parameter kn defined by Eq. (j4]). In the PCQM model {N\uu\N) = 
2 J d^rip{r)^p{r) = 1.08. 
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For the other structures of 4m condensate we have just 

(h^Ui = I f/V(v^(r)r^V^(r))'. (153) 
For the 2d2u condensate we can write 

{r^'^U = {D\dT''d\D){U\uT^u\U). 

Here \D) is the vector of state of constituent D quark. Analysis, similar to that, made for the 
4m condensates provides for the scalar- vector case X = S,Y = V 

{r^^)^^i = 2(0|Jd|0)(iV|M7oM|iV) +2 J dh-tp{r)ip{r)tp{r)-fo^{r). (154) 
Note that {N\u'you\N) = 2. For other structures 

{r^'^'Ui = J rfV(V'(r)r^^(r))'. (155) 

7.4.2 Contribution of the sea quarks 

Now all the quark operators act on the pions. The contribution is expressed in terms of pion 
expectation values of the four-quark operators 

(Il4,)p^ons = — ( ^ (^1 /^^Ur^MMF^M + T^rfF^^MF^Mlvr")) (156) 

Here a denotes the pion isotopic states, S stands for the sum of the self-energy pion loop and 
the pion exchange contribution. The pion expectation values can be expressed in terms of the 
vacuum expectation values of the four-quark operators. This was done in |125] by employing 
the current algebra technique. Using the results of |125] we find 

J^i-^^ll^xuT^uuT^u + TxdT^duT^u\iT'') = 0, (157) 

X 

and thus 

(n4,)pio„s = 0. (158) 

Hence, there is no such thing as "the contribution caused by the sea quarks only" . This is 
true for any model where the sea quarks are contained in the pions. Note that in the calculations 
of [125] the factorization assumption Eq. fll40p have been used. 

7.4.3 Contribution of the interference terms 

In the case of the scalar condensate the interference terms can be written as 

{hl)j = 2Y,{N\H,,,t\N,,7;){N,\uu\N){n\m\7f){N,^ (159) 

i 

where Hint is the interaction between the constituent quark Q and the pion. One can write 
Hint = —Lint, with Lint givcu by Eq. (11481) . Here two quark operators act on the pion, while 
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the other two operators act on a constituent quark. The latter can be the same as one in the 
matrix element of the interaction Hint or another one. Actually in the sum over the states of the 
constituent quarks Ni only those, corresponding to the nucleon are included, i.e. \Ni) = |A^). 

In the case of pseudoscalar and axial operators the interference can take place in the first 
order of the txQ interaction. This happens because the matrix elements (7r|gr^g|0) have nonzero 
values in these cases. The contribution of such "vertex interference" is 

{h^)j = (iV|if,„,|iV,7r)(iV,7r|gr^ggr^g|iV) + (iV, 7r|gr^ggr^g|iV) (iV|iJ,„i|iV, tt), (160) 

with X labelling an axial or pseudoscalar. Thus in one of the vertices interaction Hint is replaced 
by the amplitude of injection of four quarks and antiquarks. For neutral pions the latter can 
be evaluated as 

{N,Ti%T^qqT^q\N) = (7r°|gr^g|0) (iV|gr^g|iV) , 

{N\qT^qqT^q\N,n^) = {0\qT^q\7r°){N\qT^q\N), 

with analogous relations for the charged pions. 

Now we shall employ the values of the four-quark condensates, obtained in [21] in the SR 
equations. 

8 Contribution of the higher order terms 

8.1 Symmetric matter with the four-quark condensates 

We can write the general equation for the contribution of the 4q terms 

(n)4,= (X4i + xf^i^^ + x4i^4)-^p. (161) 

Note that a = — (27r)^(0|'UM|0) ~ 0.55 GeV'^ (Eq. (^Mj) is just a convenient scale for presentation 
of the results. It does not reflect the chiral properties of (11)4^. 

The coefficients Xl^ are obtained by using the complete set of the nucleon four-quark con- 
densates [3T] and the results of the previous Section. In calculations of the factor X^^ all 
contributions are numerically important due to their partial cancellations. The coefficient 
is determined mainly by the 2d2u vector-vector condensate. The factor X^^ is dominated by 
the vector-scalar 2d2u condensate, with the first term on the RHS of Eq. fll54p providing the 
largest contribution. 

The numerical values are 

Xlg = -0.11; Xfg = 0.57; = 1.27 

We took into account the nonlocal structure of the vector condensate in the factor X^^. 

Now we estimate the ratio of contributions of the 4q condensates to those of the condensates 
with d = 3, provided by Eqs. (!65|) . (!66|) . The effective values of momenta are ~ 1 GeV^. 
Putting logarithmic factors in Eqs. ( |65|) . ( |66|) equal to unity, we find that this ratio is less 
than 0.1 in the q structure of the QCD sum rules. It is about 0.13 in the P structure. The 
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ratio is about 1/4 in the scalar structure I for = 8. It becomes smaller for larger values 
of K]\f. As we shall see below, the gas approximation result for {X^)4g provided by Eq. fll62p 
overestimates the value. Thus, the values of contributions of 4g condensates are consistent with 
the assumption on the convergence of the OPE series. 

The contributions of dimension ci = 6 to the sum rules - see Eqs. (1571) . flMj) are 



2 

Ae = -8Ti^aXLp; Pg = -Svr^a X[n; Bq = -Sn^maXi^p . (163) 

^ 2m ^ 

Actually, Bq contains also the terms of the higher dimension, since it includes the nonlocality 
of the vector condensate. 

The medium induced 4g condensates on the LHS of the sum rules correspond the exchange 
by strongly correlated four-quark systems on the RHS. This may be a local two- meson exchange, 
with two mesons interaction with the nucleon at the same point - see Fig. 16. Another possible 
interpretation is the exchanges by the four-quark mesons, if there are any }126 ]. 

However, the leading contribution to the scalar structure, determined by the first term on 
the RHS of Eq. (11541) has another interpretation. Here two quarks are exchanged with vacuum, 
while two other ones are exchanged with the valence quarks of the matter. The contribution of 
this scalar- vector condensate to Xj^ is 



2{Pq) {0\dd\0)v{p) 
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This term contains the expectation value (0|(i(i|0) as a factor. On the other hand this term is 
proportional to the vector condensate, contributing, however, to the scalar Lorentz structure 
of the nucleon equation of motion Eq. ([1]). On the RHS of the SR this can be interpreted 
as the vector meson exchange with the anomalous structure of the vertex. We discussed such 
contributions in the analysis of the charge-symmetry breaking forces in Subsec. 5.2 - see Fig. 13. 

Note that if we go beyond the gas approximation, the discussed term contains rather the 
expectation value {M\dd\M) = k{p)/2 instead of the vacuum value {Q\dd\Q). As we see from 
Eq. (|92|) . |k(p)/k(0)| < 1, with the reduction of about 30 — 50% at saturation density. Hence, 
the gas approximation overestimates strongly the value of Xf^. In the next to leading order 
beyond the gas approximation we must consider the contribution, in which two pairs of quarks 
of the 4g condensate go to two different nucleons. Such terms should be considered together 
with the other three-body contributions. 

As we have seen in Sec. 4, the nucleon parameters depend on the nucleon expectation 
value kat, which is known with large uncertainties. We show density dependence of the nucleon 
parameters for kat = 8 obtained in [32] in Fig. 17. In Fig. 18 we present the dependence of 
nucleon self-energies on K^r at the saturation value of density. 

For K^r = 8 we find at p = po 
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160 MeV; m* - m = -380 MeV. (165) 



Thus, the 4g condensates subtract 90 MeV from the vector self-energy and 170 MeV from 
the effective mass m*, provided by Eq. (11051) . The potential energy U = —180 MeV looks 
discouraging. However, inclusion of the radiative corrections improves the situation. 
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8.2 Radiative corrections 



The radiative corrections to the vector condensate, calculated recently in [127J appeared to be 
rather large. Thus we have to analyze the role of the radiative corrections more carefully. 

In Subsec. 2.3 we clarified the role of these corrections in vacuum SR. Here we focus on the 
terms Up, provided by medium. The contributions of the condensates of the lowest dimension 
d = 3, provided by Eq. (jHS]) are actually written in the form 

^(0) ^ 5(0) ^ p(0) 

where the anomalous dimensions 7g = 7p = 4/9,7/ = 0- Here the upper index denotes 
that all radiative corrections are neglected. The factors L~'^^ include the sum of the terms 
{aglnq'^)"'. This is called the Leading Logarithmic Approximation (LLA). Actually, in our 
previous analysis we included the contributions of condensates with d = 3 with the radiative 
corrections treated in LLA. Note that in computation of the next to leading order terms we used 
the structure functions [68] , which reproduce their moments with proper anomalous dimensions. 
We did not include the radiative corrections to the 4q condensates, since these terms were 
obtained in framework of a model, which does not contain gluons. Also, in comparison with 
the steep dependence of the terms, containing the 4g condensates, caused by the "normal" 
high dimension, the role of anomalous dimension is relatively small. 

Now we shall include the corrections of the order beyond the LLA. We present 

i3p = 4pta; B,p = Bi%; P,p = P^lhp. (167) 

Here 



U = ^, (168) 



with 

Ti = 1 + CiOs/n, 

where the second term on the RHS is the lowest order radiative correction beyond the LLA. 
While c/ = 3/2 [51], it was found in [127] that Cg = 7/2 and cp = 15/4. Thus, employing a 
straightforward estimation, we can expect a 40 — 50% change of the nucleon parameters due to 
radiative corrections. 

We shall compare the three cases. All the radiative corrections are ignored, i.e. all tj = 1. 
Corrections are included in framework of LLA, i.e. all = 1. The third possibility is that in 
addition to LLA the lowest order «<, correction is included beyond the LLA. 

We carry out the computations for kn = 8, assuming 0^(1 GeV) = 0.37. The results of 
calculations are presented in Table 1. We see the LLA corrections subtracts 70 and 50 MeV 
from the vector self-energy and the effective mass correspondingly. Inclusion of the corrections 
beyond the LLA makes the values of Sy and m* closer to those obtained with total neglect of 
the radiative corrections. The results are illustrated by Fig. 19. 

Another problem is the dependence of numerical results on the actual value of a^. The 
latter is determined by the value of Aqcd- The authors of [127] used ^^(l GeV) = 0.47, while 
as(lGeV) = 0.55 was employed in p2]. These values corresponds to Aqcd = 0.23 GeV and 
Aqcd = 0.28 GeV. These variations of as change the nucleon self-energies by several MeV, 
affecting mostly the value of the nucleon residue [35] . 
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8.3 Asymmetric matter 

8.3.1 Inclusion of condensates of higher dimension 

The nonlocal structure of the vector condensate is included in the same way as in the case of 
the symmetric matter. 

Employing the values of the four-quark condensates obtained in the previous Section we 
find the /3 dependence of the contributions of the four-quark condensates, defined by Eq. fll6ip 

X|g = -0.11-0.21/3; Xfg = 0.57 + 0.09/3 ; = 1.27 - 0.61/3 . (169) 

Including also the nonlocal structure of the vector condensates, we solve the SR equations. 
In the actual numerical calculations we employ the PCQM value (p = 0.54 for the isotope 
asymmetric condensate of dimension d = 3 - see Eq. (IllOp . In Fig. 20, we show the values 
of the proton and neutron self-energies in neutron matter (/3 = 1), compared to thouse in 
symmetric matter. In other words, these are the proton self-energy values for the nuclear 
matter with /3 = — 1, 0, 1 [33]. 

One can see that the /3 dependence of self-energies has the same qualitative features as 
that determined by Eq. (I112p . In the matter with neutron excess /3 > the neutron vector 
self-energy is larger than the proton one. The proton effective mass is larger than the neutron 
one. In both cases the P dependence is relatively weak. 

8.3.2 Comparison with results of nuclear physics 

Our results for the difference of the effective masses of neutron and proton m*^ = m* — m* 
appears to be twice smaller than the result of |128] . but only 30% smaller than that of jl29] . 
However, a discrepancy from the relativistic Brueckner-Hartree-Fock (RHFB) calculations [129] 
is smaller. Their results for /3 = 0.2 are — = 30 MeV and m*^ = —15 MeV. We obtained 
20 MeV and —15 MeV correspondingly for these parameters. Another RHFB calculation |130] 
provided results, which are very close to ours one. Their values for /3 = 0.75 are Hy — Tiy = 
80 MeV and m*^ = —50 MeV, while we found 80 MeV and —55 MeV correspondingly. 

Note that the nuclear physics calculations provide m*^ < for /3 > 0. Our values have 
the same sign after the four-quark condensates are included. The lowest dimension solution 
expressed by Eq. (I112p provides m*^ > 0. 

We found unexpectedly satisfactory agreement for the potential energy splitting Unp = 
_ [/(p) and even for the /3 dependence of the average binding energy per nucleon e{p,P). 
Various approaches |13H 1132] provided Unp ~ 60 MeV, while our result is Unp = 40 MeV. The 
"symmetry energy" e^^"^ = l/2de{po, f3)/df3 is 29 MeV in our approach. All earlier |128j - [T^ 
and nowadays |133j calculations provide e^^"^ ~ 30 MeV. We did not expect a good agreement 
here, since these parameters contain the large terms, which cancel each other to large extent. 

8.4 Many-body interactions 

The structure and the role of the three-body (many-body) nuclear interactions is much studied 
nowadays - see, e.g. [HI I134[ 1135] . In the SR approach such interactions emerge in a natural 
way. 
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There are two types of many-body interactions in our approach. One of them is connected 
with the four-quark condensate. Another one manifests itself in the baryon-hole excitations, 
taken into account by the in-medium pion propagator. 

In Sec. 6 the contribution of four-quark condensates was included in the gas approximation. 
Our probe nucleon exchanged both pairs of quarks with the same nucleon of the matter. Beyond 
the gas approximation two quark pairs can be exchanged with different nucleons of the matter. 
This corresponds to the three-body interactions. 

Such terms are shown in Fig. 21. Their contribution, e.g. to the LHS of Eq. (!57|) for the q 
structure C^i^M"^ ,W^) can be estimated as Att'^k^p^. At the saturation value of density this 
can change the value of £^(M^, W^) by about 10%. This will lead to relative changes of the 
same order of the nucleon self-energies. At larger values of density the role of these terms 
increases. 

As we discussed in Sec. 3, the contributions, involving n > 2 nucleons of the matter have 
the branching points in variable Sniq^) = {nP + g)^. Generally speaking, this means that we 
need a more complicated model of the spectrum. However, we see that in the OPE series the 
three body terms do not contain branching points in g^. In framework of our model for the 
higher lying states they contribute only to the pole term on the RHS of the sum rules. 

Adding gluon interactions with the matter, i.e. including the four-quark-gluon condensates, 
we obtain the terms, involving larger number of nucleons of the matter, corresponding to many- 
nucleon forces. They contribute to the higher order OPE terms. Due to the small value of the 
in-medium gluon condensate we expect such contributions to be small. 

Another type of many-body interactions manifests itself in the nonlinear contribution to 
the scalar condensate S{p). As we saw in Subsection 4.5, it is determined by the pion cloud. 
The general form of the contribution is 

- Tf(k)Df«y\k)Gf(p-k)Tf(k)){,r\qq\,r), (170) 

where summation over spin and isospin variables is assumed. Here is the propagator 
of intermediate baryon in nuclear matter, D^'^\k) is the pion propagator, renormalized by 
the particle-hole excitations of medium. The vertices of Nt[B interaction in nuclear matter 
are denoted as P^. Integration over momenta of the nucleons of the matter is limited by 
condition p < pp. Upper index (0) denotes the vacuum functions. The second term on the 
RHS of Eq. (11701) subtracts the terms, which are already included in the physical nucleon. For 
j^im) ^ ^(0). = the RHS of Eq. ([ITOD turns to that of Eq. ([IISD. The contribution is 
illustrated by Fig. 22. 

The pion propagator D^'^\k) includes the multi-nucleon effects. Thus, Eq. fll70p includes 
the many-body interactions. However, the rigorous calculation of its RHS is a part of a self- 
consistent problem. 



47 



9 Self-consistent scenario 



We have seen that the shape of the density dependence of the quark condensates ^(p) is very 
important for the hadronic physics. It is also important for description of the matter as a whole, 
since characterizes a degree of restoration of the chiral symmetry with the growing density. On 
the other hand, we saw that it determines the main density dependence of the effective mass 
m*(p) ~ m*(K(p)). 

The linear part of the density behavior of k{p) is related to observables and can be calculated 
in a model-independent way. However, the nonlinear contribution, expressed by Eq. (I170p 
depends on several hadron parameters. If the intermediate baryon in Eq. fll70p is a nucleon, 
these are the nucleon effective mass m* and the nucleon pion coupling constant QnNNip)- The 
latter can be expressed in terms of the pion decay constant fjrip) and the nucleon isovector 
axial coupling qaIp) via the Goldberger-Treiman relation [136\ 

g-KNN/'^rn = QA/fn- 

The density dependence of qa and /,r can also be obtained in the SR approach. One should add 
similar equations for the contribution of delta- isobars to the RHS of Eq. fll70p . Thus we come 
to self-consistent set of equations which can be solved within the QCD sum rules approach. 

m*=m*{K); = /^(k); gA = 9AiK,); n = n{m\ f^/ qa). (171) 

Of course, the SR approach should be combined with some other ones. We have seen already 
that the calculation of the four-quark condensate required a quark model for the nucleon. Also, 
the RHS of Eq. (11701) contains the in-medium pion propagator. It satisfies the Dyson equation 

where the pion polarization operator H^r describing the baryon-hole excitations of the matter. 
The operator H^^ depends strongly on the hadron correlations at the distances r > 
corresponding to small momenta. They can be included by employing of the Finite Fermi 
System Theory (FFST) |137j . [138] . Thus, in order to obtain the many-body effects it is 
reasonable to combine the SR approach with the FFST. 

10 Summary 

We demonstrated that the QCD sum rules provide a consistent formalism for solving various 
problems of nuclear physics. 

We expressed the nucleon parameters in nuclear matter in terms of the in-medium QCD 
condensates. In symmetric matter the leading contributions to the vector and scalar self- 
energies are expressed in terms of vector and scalar condensates. The vector condensate v{p) 
can be calculated easily. In the gas approximation, corresponding to inclusion only of the 
two-body interactions, the scalar condensate k[p) is related to the observable pion-nucleon 
sigma-term. 

In other words, exchange by the strongly correlated quark systems (mesons) is expressed in 
terms of exchange by the system of weakly interacting quarks with the same quantum numbers. 
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Thus the nucleoli self-energies are obtained without employing a controversial conception of 
interaction of the point-like nucleons. 

Inclusion of the condensates of the lowest dimension confirmed that the nucleon in nuclear 
matter can be treated as moving in superposition of strong vector and scalar fields of the order 
of about 300 MeV. The vector field is positive, while the scalar field is negative, and a strong 
cancellation takes place. 

This result is not altered by inclusion of the four-quark condensates. Thus we can expect the 
convergence of the series with successive inclusion of the condensates of higher dimension. Note, 
however, that the calculation of the in-medium four-quark condensates requires employing a 
quark model of nucleon. 

We used the approach for calculation of other nucleon characteristics in nuclear matter. We 
calculated the nucleon self-energies in asymmetric nuclear matter and found their dependence 
on the difference between the densities of the neutrons and protons. We calculated the in- 
medium quenching of isovector axial coupling constant. We found the neutron-proton mass 
difference, providing at least qualitative explanation of the Nolen-Schiffer anomaly. 

Also, we demonstrated that our approach enables to investigate internal structure of the nu- 
cleon, providing description of some features of the EMC effect. Such problems are unaccessible 
for nuclear physics. 

Inclusion of the nonlinear density dependence of the scalar condensate K,{p) corresponds to 
taking into account the many-body interactions. We show that the nonlinear terms in K,{p) 
are provided mainly by the pion cloud. A simple inclusion of the contributions beyond the 
gas approximation signals on a possible saturation mechanism, which is due to the many-body 
effects. Thus it differs from that in the Walecka model. 

A more rigorous treatment of the nonlinear contributions to /t(p) requires solution of the 
self-consistent problem, in which the coupling constants Qt^nn and g.,rNA as well as the effective 
mass of the A isobar should be also determined as functions of k, while k should be expressed 
in terms of these parameters. This would provide the dependence of k{p). This is an important 
parameter of the matter itself, since it shows the degree of restoration of the chiral symmetry. 

This analysis requires employing of the pion propagator, which is renormalized by the 
particle-hole excitations of the matter. These excitations can be included by using the Finite 
Fermi System Theory, formulated and developed by A. B. Migdal and his colleagues several 
decades ago. 

The work on the project is in progress. 

11 Epilogue 

Many years ago the authors of this paper participated in a very interesting discussion on the 
perspectives of nuclear physics, which took place at the traditional Petersburg (Leningrad at 
that time) Winter School of Physics. In his lecture A. B. Migdal formulated two questions, he 
wished to be answered in the future: 

• How do the QCD condensates change in nuclear medium ? 
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• What is the connection between the in-medium QCD condensates and hadron parame- 
ters? 

We hope that the present paper makes the first steps to answer these questions. 
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Table 1: Nucleoli parameters at the saturation value of nucleon density. Line 1: all radiative 
corrections are neglected. Line 2: radiative corrections are included in LLA. Line 3: corrections 
~ CKs are included beyond the LLA (BLLA). 







m* — m 


\2 


' ' m 




(MeV) 


(MeV) 


(GeV6) 


(GeV2) 


No corrections 


229 


-329 


1.10 


1.72 


LLA 


160 


-380 


1.10 


1.77 


BLLA 


271 


-300 


1.41 


1.75 



Table 2: Dependence of the nucleon parameters on the values of Aqcd and ag at the phe- 
nomenological saturation value of nucleon density with radiative corrections included beyond 
the leading logarithmic approximation (BLLA). The results are presented for Aqcd = 0.23 GeV, 
Q;^(lGeV^) = 0.47, and Aqcd = 0.28 GeV, ^^(IGeV^) = 0.55. The corresponding values of 
vacuum parameters arc given in brackets. 
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(1 GeV2) 
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(GeV^) 


0.47 


269 
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Figure Captions 

Fig. 1. Main contributions to the LHS of the vacuum sum rules - Eq. (19). Helix line stands 
for the correlator, solid lines denote the quarks, dotted lines stand for the quarks composing 
the scalar condensate. Figs, a, b and c correspond to the contributions Aq, B3 and Aq. 

Fig. 2. A typical radiative correction. The dashed line denotes the gluon. 

Fig. 3. Contributions to the singularities of the correlator in s channel (a) and in u channel 

(b). 

Fig. 4. Self-energy insertions to the nucleon pole contributions: the mean field approxima- 
tion (a), the self-energy in the direct channel (6) and the exchange contribution (c). 

Fig. 5. One of the contributions, containing a higher lying singularity in at fixed s. Helix 
line stands for the correlator. Dashed lines stand for the mesons, solid line is for the nucleon 
and thick the solid line denotes the nuclear matter. 

Fig. 6. A singularity in u channel, containing an intermediate state with the baryon number 
B = 0. The meaning of the lines is the same as in Fig. 5. 

Fig. 7. Contributions of the condensates of lowest dimension to the left-hand side of the sum 
rules in nuclear matter. The dotted lines in (a) stand for the quarks of the scalar condensate. 
The dashed-dotted lines in (6) denote the quarks of the vector quark condensate. The meaning 
of the other lines is the same as in previous figures. 

Fig. 8. Contributions of the condensates of lowest dimension to the right-hand side of the 
sum rules in nuclear matter. The horizontal line denotes the nucleon, thick lines show the 
matter. The wavy and dashed lines stand for vector and scalar mesons. 

Fig. 9. The dependence of the functions fp and // defined by Eq. ( I5T]) on the Borel 
mass M^. 

Fig. 10. Density dependence of the nucleon parameters for = 8: the vector self-energy 
and effective mass (a); the residue and the threshold (6), related to their vacuum values. 
The X axis corresponds to the density, related to its saturation value. Solid line is for solution 
of Eq.(57) with condensates of dimensions d = 3 and d = 4 taken into account. Dashed lines 
show the approximate solution, corresponding to Eqs. (89)-(91). 

Fig. 11. Dependence of the nucleon parameters on the value of kn at saturation density. 
Solid line is for the effective mass m*, dotted line shows the vector self-energy Ey, dashed and 
dash-dotted lines are for the and correspondingly. All parameters are related to their 
vacuum values (the vector self-energy is related to the vacuum value of the nucleon mass) given 
by Eq.(^. 

Fig. 12. Simplest contribution to the nonlinear scalar condensate. Solid lines labeled as 1 
and 2 are the nucleons of the matter. Dashed line denotes a vr meson. 

Fig. 13. Anomalous Lorentz structure of nucleon interactions with matter. Notation are 
the same as in Fig. 8. 

Fig. 14. Second order interaction of the hard photon (dashed lines) with the correlator. The 
dark blobs denote interaction with the matter. Other notations are the same as in Fig. 5. 

Fig. 15. The in-medium changes of the d quark distribution (dashed curve) and of the u 
quark distribution (dot-dashed curve) of the fraction x of the momentum of the target nucleon. 
The solid curve presents the function R — 1 with the ratio R, defined by Eq. (126). 

Fig. 16. Interactions on the RHS of the sum rules, corresponding to inclusion of the four- 
quark condensates. The wavy lines denote mesons. 



Fig. 17. Density dependence of the nucleon parameters for — 8 with the condensates of 
higher dimensions included: the vector self-energy and effective mass (a); the residue A^^ and 
the threshold (b), related to their vacuum values. The x axis corresponds to the density, 
related to its saturation value. 

Fig. 18. Dependence of the nucleon parameters on the value of kjv at saturation density. 
Condensates of higher dimension are included. The meaning of the curves is the same as that 
for Fig. 11. 

Fig. 19. Density dependence of the vector self energy and of the scalar self energy 
m* — m. The nuclear matter density p is related to its saturation value po- Dotted lines: 
all radiation corrections are neglected. Dashed lines: radiative corrections are included in the 
Leading Logarithm Approximation (LLA). Solid lines: Corrections of the order ccg are included 
perturbatively beyond the LLA. 

Fig. 20. The density dependence of the vector self-energy Sy (a) and of the effective mass 
m* (6) in isospin asymmetric matter. The solid line is for symmetric matter. The dashed and 
dotted lines are for the proton and neutron characteristics in neutron matter (/3 = 1). 

Fig. 21. A contribution of the three-body forces to the LHS of the sum rules. Two pairs 
of quarks from the four-quark condensate interact with two different nucleons of the matter, 
shown by the dotted circles. 

Fig. 22. Interaction of the operator qq (the dark blob) with pion field. The solid line denotes 
the nucleon of the matter; the wavy line stands for the pion. The bold wavy hne denotes the 
pion propagator renormalized by the baryon-hole excitations. 
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